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THE COMPRESSIBLE VISCOUS SURFACE-INTERNAL WAVE PROBLEM: 

LOCAL WELL-POSEDNESS 

JUHI JANG, IAN TICE, AND YANJIN WANG 


Abstract. This paper concerns the dynamics of two layers of compressible, barotropic, viscous 
fluid lying atop one another. The lower fluid is bounded below by a rigid bottom, and the upper 
fluid is bounded above by a trivial fluid of constant pressure. This is a free boundary problem: 
the interfaces between the fluids and above the upper fluid are free to move. The fluids are 
acted on by gravity in the bulk, and at the free interfaces we consider both the case of surface 
tension and the case of no surface forces. We prove that the problem is locally well-posed. 
Our method relies on energy methods in Sobolev spaces for a collection of related linear and 
nonlinear problems. 


1. Introduction 

1.1. Formulation in Eulerian coordinates. We consider two distinct, immiscible, viscous, 
compressible, barotropic fluids evolving in a moving domain Q(t) = H+(t) U fl_(t) for time 
t > 0. One fluid (+), called the “upper fluid,” fills the upper domain 

M+(t) = {y € T 2 x M | r]-(yi,y 2 ,t ) < y 3 < £ + rj + (y 1 ,y 2 ,t)}, (1.1) 

and the other fluid (—), called the “lower fluid,” fills the lower domain 

O-(t) = {y € T 2 x R | -b < y 3 < y-(yi, y 2 , t)}. (1.2) 

Here we assume the domains are horizontal periodic by setting T 2 = ( 2-kL\T ) x ( 2itL 2 T ) for 
T = M/Z the usual 1-torus and L\,L 2 > 0 the periodicity lengths. We assume that £,b > 0 
are two fixed and given constants, but the two surface functions y± are free and unknown. The 
surface r + (f) = {y 3 = £ + r] + (yi, y 2 ,t)} is the moving upper boundary of H + (t) where the upper 
fluid is in contact with the atmosphere, T_(t) = {y 3 = ri_(yi,y 2 ,t)} is the moving internal 
interface between the two fluids, and E& = {2/3 = —6} is the fixed lower boundary of fl_(t). 

The two fluids are described by their density and velocity functions, which are given for each 
t > 0 by p±(-,t) : £l±(t) -A M + and u±(-,t) : fi±(t) —> M 3 , respectively. In each fluid the 
pressure is a function of density: P± = P±(p±) > 0, and the pressure function is assumed to be 
smooth, positive, and strictly increasing. For a vector function u € M 3 we define the symmetric 
gradient by (Dii)^ = diUj + djUi for i,j = 1,2,3; its deviatoric (trace-free) part is then 

n 2 

D°tt = Ou — - div ul, (1.3) 

O 

where I is the 3x3 identity matrix. The viscous stress tensor in each fluid is then given by 

§±(h±) := yU-i-B 0 ^ + //()_ div u±I, (1.4) 

where y± is the shear viscosity and p'± is the bulk viscosity; we assume these satisfy the usual 
physical conditions 

y± > 0, y'± > 0. (1.5) 

The tensor P±(p±)I — §±(h±) is known as the stress tensor. The divergence of a symmetric 
tensor M is defined to be the vector with components (divM)j = djWlij. Note then that 

di v(P±(p±)I ~ S±(h±)) = VP±(p±) - p±Au± - (J^ + p'±J V div u±. (1.6) 
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For each t > 0 we require that (u±, p±,rj±) satisfy the following equations: 


'dtp± + div(p±u±) = 0 in £l±(t) 

p±{d t u± + u± ■ Vu±) + VP±(p±) - divS±(it±) = -gp±e 3 in Q±(t) 
<9tp± = «3,± - «i,±9 yi ??± - u 2) ±dy 2 r]± on r±(t) 

< (P+(p + )J - § + (£t+))n + = Pat m n + - a + n + n + on T + (t) 


(P+(p+)I — § + (£t_|_))n_ = (P_(p_)7 — §_(u_))n_ + cr_'H_n_ on r_(f) 
£t_l_ = U- on r_(t) 

= 0 on Sft. 


(1.7) 


In the equations —gp±e 3 is the gravitational force with the constant g > 0 the acceleration of 
gravity and e% the vertical unit vector. The constant p a tm > 0 is the atmospheric pressure, 
and we take a± > 0 to be the constant coefficients of surface tension. In this paper, we let 
V* denote the horizontal gradient, div* denote the horizontal divergence and A* denote the 
horizontal Laplace operator. Then the upward-pointing unit normal of T±(t), n±, is given by 


n± 


(-V*q±. !) 

y/l+WJhP’ 


( 1 . 8 ) 


and 'H±, twice the mean curvature of the surface T±(t), is given by the formula 


U± 


div* 


( V*7?± \ 

Ul + I V *M 2 / ' 


(1.9) 


The third equation in (11.711 is called the kinematic boundary condition since it implies that the 
free surfaces are advected with the fluids. The boundary equations in (11,71) involving the stress 
tensor are called the dynamic boundary conditions. Notice that on T_(f), the continuity of 
velocity, = u_, means that it is the common value of u± that advects the interface. For a 
more physical description of the equations (11.71) and the boundary conditions in (11.71) . we refer 
to [ 32 ]- 

To complete the statement of the problem, we must specify the initial conditions. We suppose 
that the initial surfaces r(0) are given by the graphs of the functions ^±(0), which yield the 
open sets H±(0) on which we specify the initial data for the density, p±(0) : H±(0) —>• M + , and 
the velocity, fi±(0) : fi±(0) —> M 3 . We will assume that i + 77 +(0) > r/_(0) > — b on T 2 , which 
means that at the initial time the boundaries do not intersect with each other. 


1.2. Equilibria. Here we seek a steady-state equilibrium solution to m with u± = 
and the equilibrium domains given by 

= {y € T 2 x R | 0 < 2/3 < 1} and = {y € T 2 x R | — b < 7/3 < 0}. 


Then (11.71) reduces to an ODE for the equilibrium densities p± = p±(yz)'- 


d(P+(p +)) 

dy 3 

d{P-(p-)) 

dm 

d D +{P+{^)) = Patm: 

P+(p+(0))=P_(p_(0)). 


= -9 P+, 


= -gp-, 


for 2/3 € ( 0 , £), 
for 2/3 € (-6,0), 


0, V± = 0, 

( 1 . 10 ) 


( 1 . 11 ) 


The system (11.111) admits a solution p± > 0 if and only if we assume that the equilibrium 
heights b,£ > 0, the pressure laws P±, and the atmospheric pressure p a tm satisfy a collection of 
admissibility conditions. These are enumerated in detail in our companion paper M- For the 
sake of brevity we will not mention these here, but we will assume they are satisfied so that an 
equilibrium exists. The resulting function p is strictly positive and smooth when restricted to 
to [—6,0] and [0, t\. 

We give special names to the equilibrium density at the fluid interfaces: 

Pi = P+{£), P + = /M0), P~ = M0). 


( 1 . 12 ) 
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Notice in particular that the equilibrium density can jump across the internal interface. The 
jump in the equilibrium density, which we denote by 

M ■= P+(0) ~ P-{0) = P + ~ P~, (1-13) 

is of fundamental importance in the the analysis of solutions to (11.71) near equilibrium. Indeed, 
if IpI > 0 then the upper fluid is heavier than the lower fluid along the equilibrium interface, 
and the fluid is susceptible to the well-known Rayleigh-Taylor gravitational instability. This 
is not particularly important for the local theory developed in this paper, but of fundamental 
importance in the stability theory. 

In studying perturbations of the equilibrium density it will be useful throughout the paper 
to employ the enthalpy functions. These are defined in terms of the pressure laws P± and the 
equilibrium density values via 

r~ P' ( r ) r z P' (r) 

h + (z ) = / ——-dr and h-(z ) = / —— -dr. (1-14) 

J pi r Jp- r 

1.3. Reformulation in flattened coordinates. The movement of the free surfaces T-j-(f) 
and the subsequent change of the domains create numerous mathematical difficulties. 

To circumvent these, we will switch to a coordinate system in which the boundaries and the 
domains stay fixed in time. In order to be consistent with our study of the nonlinear stability of 
the equilibrium state in m, we will use the equilibrium domain. We will not use a Lagrangian 
coordinate transformation, but rather utilize a special flattening coordinate transformation 
motivated by Beale [3] . 

To this end, we define the fixed domain 

fl = 11+ U with := {0 < < £} and := {—6 < x 3 < 0}, (1-15) 

for which we have written the coordinates as x £ Q. We shall write £ + := {£3 = £} for the 
upper boundary, £_ := { 2:3 = 0 } for the internal interface and £& := { 2:3 = — 6 } for the lower 
boundary. Throughout the paper we will write £ = £_|_ U £_. We think of r/± as a function on 
£± according to r] + : (T 2 x {£}) x R + — > R and : (T 2 x {0}) x R + R, respectively. We 
will transform the free boundary problem in fl(t) to one in the fixed domain fl by using the 
unknown free surface functions T]±. For this we define 

f) + := V+r] + = Poisson extension of r/ + into T 2 x { x 3 < £} (1-16) 

and 

77 _ := V-T = specialized Poisson extension of into T 2 x R, (1-17) 

where V± are defined by (1C.3D and (IC. 8 D . The Poisson extensions fj± allow us to flatten the 
coordinate domains via the following special coordinate transformation: 

9 x i-A ( x 1 ,x 2 ,x 3 + bifj + + 6 2 i?-) := = ( 2 / 1 , 2 / 2 , 2 / 3 ) € (1.18) 

where we have chosen b\ = 61 ( 2 : 3 ), 6 2 = 62 ( 2 : 3 ) to be two smooth functions in R that satisfy 

6 i( 0 ) = 6 i(- 6 ) = 0 , 61 (f) = 1 and b 2 {£) = 6 2 (—6) = 0,6 2 (0) = 1. (1-19) 

Note that 0(£ + ,t) = r+(t), 0(£_,t) = r~(t) and 0(-,t) |s 6 = Id |s 6 . 

Note that if r/ is sufficiently small (in an appropriate Sobolev space), then the mapping 0 
is a diffeomorphism. This allows us to transform the problem m to one in the fixed spatial 
domain for each t > 0. In order to write down the equations in the new coordinate system, 
we compute 

/ 1 0 0 \ / 1 0 -AK \ 

V0 = I 0 1 0 and A := (V0 _1 ) T = I 0 1 -BK . (1.20) 

\ABjJ \ 0 0 K J 

Here the components in the matrix are 

(1.21) 


A = <9i 6, B = d 2 0, J = 1 + d 3 0 , K = J 
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where we have written 

9 :=bifj+ + b 2 fj-. ( 1 - 22 ) 

Notice that J = detV0 is the Jacobian of the coordinate transformation. It is straightforward 

to check that, because of how we have defined p- and 0, the matrix A is regular across the 

interface £_. 

We now define the density p± and the velocity u± on Q± by the compositions p±(x,t ) = 
p±(@±(x,t),t) and u±(x,t) = u±(Q±(x,t),t)- Since the domains and the boundaries £± 
are now fixed, we henceforth consolidate notation by writing / to refer to f± except when 
necessary to distinguish the two; when we write an equation for / we assume that the equation 
holds with the subscripts added on the domains or To write the jump conditions on 
for a quantity f = f±, we define the interfacial jump as 

Ifl ;= f+\{x 3 =0} — /-l{a: 3 = 0 }- (1-23) 

Then in the new coordinates, the PDE (11.71) becomes the following system for ( u,p,rj ): 


dtp — Kdtddsp + d\\ A (pu) = 0 in 17 

p(d t u - Kd t 9d 3 u + u • X A u) + X A P(p) - div^S^(u) = -gpe 3 in Q 

dtp = u ■ M on E 

< ( P(p)I - (u)W = PatmN - V+UN on S + 

\P{p)I — §a( w )1 -A/" = oAKM on S_ 

[uj = 0 on E_ 

^U- = 0 on E{,. 


(1.24) 


Here we have written the differential operators V_ 4 , div_ 4 , and A _4 with their actions given by 
(V A f)i ■= Ajdjf, div_ 4 X := AijdjXi, and A A f := div A V A f (1.25) 

for appropriate / and X. We have also written 

AT := {-d 1 p,-d 2 pA) (1-26) 

for the non-unit normal to S(t), and we have written 


n Z, 

(ID ) A u)ij = AikdkUj + AjkdkUi, D A u = D A u - - div^ ul , 

and §a,±( w ) := p±^> A u + p± div^ ul. (1-27) 

Note that if we extend divyi to act on symmetric tensors in the natural way, then div _4 S A u = 
pA A u + (p/3 + p')X A div _4 u. Recall that A is determined by p through (11.201) . This means 
that all of the differential operators in (ll.24|) are connected to p, and hence to the geometry of 
the free surfaces. 


Remark 1.1. The equilibrium state given by (11.111) corresponds to the static solution (u,p,p) = 

(o,p,o) of (HMD. 

1.4. Perturbation equations. We will now rephrase the PDE (11.241) in a perturbation for¬ 
mulation around the steady state solution (0, p, 0). We define a special density perturbation 
by 

q = p-p-d 3 p0. (1.28) 

In order to deal with the pressure term P(p) = P(p + q + d 3 p0) we introduce the Taylor 
expansion: by (11.111) we have 

P(p + q + d 3 pd) = P(p) + P'(p)(q + d 3 p0) +U = P(p) + P'(p)q - gp9 + K, (1.29) 
where the remainder term is defined via 

rp+q+dspG 

1Z= (p + q + d 3 p6 — z)P"(z) dz. (1.30) 

Jp 
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The advantage of defining the perturbation in this manner is seen in the following cancellation: 

AijdjP(p) + gpdi 3 = AijdjP(p) + gpAijdjQ 3 

= Aijdj{P(p) + P'{p)q - gpO + TZ) + g(p + q + d 3 p0)Aijdj(x 3 + 6) 

= Aijdj^P'(p)q) - gAijdj(p6 ) + AijdjTZ + gpAijdjO + g(q + d 3 p8)A i3 + g{q + d 3 pd)AijdjO 
= AijdjiP' (p)q) + AijdjTZ + gqA i3 + g{q + d 3 p8)A ij d j 8 
= pAijdj(h'(p)q ) + AijdjTZ + g(q + d 3 p0)Aijdj0, 

(1.31) 


where we have used (11.111) and (|1.14l) . Recalling also (11.131) . (11.191) . and (11.221) . we find that 


- gp+9 = ~Pi9P+ on £+, and l~gp8j = - [p] gp- on £_. (1.32) 


The equations (11.241) become the following system when perturbed around the equilibrium 
(0, p, 0): 


dtq + div 4 ((p + q + d 3 pd)u) — d\pK0dt8 — Kdf0d 3 q = 0 in fl 

(,P + q + d 3 pd)d t u + {p + q + d 3 pd)(—Kd t 8d 3 u + u ■ V^rt) + pS7 a (h'(p)q) 

— div .4 = —VAJZ — g(q + d 3 p0)\Z^0 in fl 

dtp = u ■ AT on E 

( P'(p)ql — §a(u))AT = pigpM — (J+TiM — 7Z+AT on £ + 

\P'{p)qI - §a(«) 1 AC = [p] gpM + vPHN - pZj M on 

[it] =0 on 

U- = 0 on Sb- 


(1.33) 


Remark 1.2. The introduction of the special density perturbation q given by (J1.28D and the 
subsequent perturbation equations of form (11.331) is crucial for our study of the nonlinear stability 
in [H ]; it is not essential for the local theory developed in this paper. Indeed, we could consider 
p or p — p directly. We choose here to consider q in order to be consistent with the study in M- 


2. Main results and discussion 

2.1. Previous work. Free boundary problems in fluid mechanics have attracted much interest 
in the mathematical community. A thorough survey of the literature would prove impossible 
here, so we will primarily mention the work most relevant to our present setting, namely that 
related to layers of viscous fluid. We refer to the review of Shibata and Shimizu m for a more 
proper survey of the literature. 

The dynamics of a single layer of viscous incompressible fluid lying above a rigid bottom, 
i.e. the incompressible viscous surface wave problem, have attracted the attention of many 
mathematicians since the pioneering work of Beale [3]. For the case without surface tension, 
Beale [3] proved the local well-posedness in the Sobolev spaces. Hataya m obtained the global 
existence of small, horizontally periodic solutions with an algebraic decay rate in time. Guo 
and Tice mmm developed a two-tier energy method to prove global well-posedness and 
decay of this problem. They proved that if the free boundary is horizontally infinite, then the 
solution decays to equilibrium at an algebraic rate; on the other hand, if the free boundary is 
horizontally periodic, then the solution decays at an almost exponential rate. The proofs were 
subsequently refined by the work of Wu [33]. For the case with surface tension, Beale [1] proved 
global well-posedness of the problem, while Allain [1] obtained a local existence theorem in two 
dimension using a different method. Bae [2] showed the global solvability in Sobolev spaces via 
energy methods. Beale and Nishida [5] showed that the solution obtained in [3] decays in time 
with an optimal algebraic decay rate. Nishida, Teramoto and Yoshihara [18] showed the global 
existence of periodic solutions with an exponential decay rate in the case of a domain with a 
flat fixed lower boundary. Tani m and Tani and Tanaka |28| also discussed the solvability 
of the problem with or without surface tension by using methods developed by Solonnikov in 
[22 , 23, 123]. Tan and Wang [2S] studied the vanishing surface tension limit of the problem. 



















6 


JUHI JANG, IAN TICE, AND YANJIN WANG 


There are fewer results on two-phase incompressible problems, i.e. the incompressible viscous 
surface-internal wave or internal wave problems. Hataya [12j proved an existence result for a 
periodic free interface problem with surface tension, perturbed around Couette flow; he showed 
the local existence of small solution for any physical constants, and the existence of exponentially 
decaying small solution if the viscosities of the two fluids are sufficiently large and their difference 
is small. Priiss and Simonett [T9] proved the local well-posedness of a free interface problem 
with surface tension in which two layers of viscous fluids fill the whole space and are separated 
by a horizontal interface. For two horizontal fluids of finite depth with surface tension, Xu and 
Zhang [33] proved the local solvability for general data and global solvability for data near the 
equilibrium state using Tani and Tanaka’s method. Wang and Tice [30] and Wang, Tice and 
Kim m adapted the two-tier energy methods of mmm to develop the nonlinear Rayleigh- 
Taylor instability theory for the problem, proving the existence of a sharp stability criterion 
given in terms of the surface tension coefficient, gravity, periodicity lengths, and [pj. 

The free boundary problems corresponding to a single horizontally periodic layer of com¬ 
pressible viscous fluid with surface tension have been studied by several authors. Jin m and 
Jin-Padula m produced global-in-time solutions using Lagrangian coordinates, and Tanaka 
and Tani [26] produced global solutions with temperature dependence. However, to the best 
of our knowledge, even the local existence problem for two layers of compressible viscous fluids 
remains unsolved. 

The two-layer problem is important because it allows for the development of the classical 
Rayleigh-Taylor instability [20l [29], at least when the equilibrium has a heavier fluid on top 
and a lighter one below and there is a downward gravitational force. In our companion paper 
m we identify a stability criterion and prove the existence of global solutions that decay to 
equilibrium. In our companion paper we show that the stability criterion is sharp, as 
in the incompressible case 122 ED, and that the Rayleigh-Taylor instability persists at the 
nonlinear level (the linear analysis was developed by Guo and Tice in [3]). The Rayleigh-Taylor 
instability is a long time phenomenon; for the local-in-time theory developed in this paper it 
plays no essential role. 


2.2. Local existence. In our companion paper m we deal with questions of global existence 
and asymptotic stability of solutions to the perturbed system (11.3311 . The analysis there is 
carried out in a high-regularity functional framework that is indexed by an integer N > 3 
related to the decay properties of solutions. Consequently, we must first guarantee the local 
existence of solutions in this framework for every N >3. 

Our main result guarantees the existence of high-regularity solutions to (11.331) under smallness 
conditions on the initial free surface as well as the temporal interval. Notice, though, that there 
are no smallness conditions placed on the initial velocity or density. We refer to Appendix [A] 
for the definitions of the terms £q , 6, T> , £ a , T> a , and £ (which all depend on N ) appearing in 
the statement of the theorem. 

Theorem 2.1. Let N >3 be an integer. Assume that either cr± > 0 or else cr± = 0. Suppose 
that (uo,qo,r]o) satisfy £q < oo in addition to the compatibility conditions (IB. 41 ) . and that 

P* < Po ■= P + qo + d 3 pdo < p* (2.1) 

for two constants 0 < p*, p* < oo. There exists a universal constant 5^ > 0 and a Ti oc = T) oc (£q ) 
such that if 

£[%] < y and 0 < T < T ioc , (2.2) 

then there exists a triple (u, q, p) defined on the temporal interval [0, T) satisfying the following 
three properties. First, (u,q,p) achieve the initial data at t = 0. Second, the triple uniquely 
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solve (11.331) . Third, the triple obey the estimates 


sup (£[u(t)} + £[q(t)} +£ a [r)(t)] + \\v(t )\\ 4N+1/2 

0 <t<T v 


+ 


V[u{t)\ + pJd^ N+ 1 u(t) 


(oH 1 )- 


+ V[q(t)}+V"[r,(t)} )dt<P(Sg) (2.3) 


for a universal positive polynomial P with P( 0) = 0. Also, 


and 


1 3 

-p* < P( x -> := P( x ) + q( x i t) + 83 p(x)9(x, t) < -p* for all x € 12 and t € [0, T] 


m(T)<6 v . 


(2.4) 


(2.5) 


Moreover, the mapping @(-,i) defined by ([1.181) is a C 4N 2 diffeomorphism for each t € [0, T\. 

Remark 2.2. The diffeomorphism condition in Theorem \2.1\ guarantees that the solution to 
(11.331) also gives rise to a solution to (11.71) in the original coordinate system. 


Remark 2.3. The temporal existence interval of Theorem \2. 1\ with u± > 0 is independent of cr±. 
As such, we can use a standard limiting argument to produce a solution to (11.331) when <t+ > 0 
and a- = 0 or a + = 0 and (j_ > 0 (the case a± = 0 is already covered by the theorem). While 
this produces a solution, it does not yield a satisfactory local well-posedness theory. Indeed, if 
we are given data (uo, go, t?o) satisfying the compatibility conditions (IB.41) with, say cr+ > 0 and 
(7 — = 0, we would need to produce a sequence of approximate data (u^qQ,^) satisfying ()B.4D 
with cr_j_ > 0 that converge to (uq, qo, 770) as ct_ — >■ 0. While this construction might be possible, 
the compatibility conditions El) are sufficiently complicated that we have not pursued it in this 
paper. 

It is noteworthy, though, that if we send both a± —>• 0, then the same sort of limiting argument 
can provide a sort of continuity result, connecting our result with a± > 0 to our result with 
a± = 0. See [Hj for more details. 


Theorem no is proved in a somewhat more general form (the smallness conditions on the 
data are more general) in Theorems 19.71 and 110.11 The former handles the case a± > 0, while 
the latter handles the case cr± = 0. The proof of Theorem 19.71 is more complicated than that of 
Theorem llU.il because of the regularity requirements for dealing with the surface tension terms 
in (11.331) . Therefore, most of the paper is devoted to the proof of Theorem 19.71 We sketch below 
some of the main ideas of the proof. 

Picard iteration 

We construct solutions to (11.33|) with a± > 0 by way of a Picard iteration scheme, which is 
developed in Section [9l This scheme is built on two sub-problems: the transport problem (18.11) 
for q, and the two-phase free boundary Lame problem (17.11) for ( u,rj ). The sequence of Picard 
iterates {(u n , q n , rj n )}ff =0 is constructed in Theorem 19.51 under some smallness assumptions on 
r/o and time. The theorem also provides estimates of various high-regularity norms in terms of 
polynomials of . 

The uniform bounds are sufficient for extracting weak and weak-* limits from the sequence of 
Picard iterates. Space-time compactness results then provide strong convergence results. How¬ 
ever, the nature of our iteration scheme does not allow us to immediately deduce from the strong 
subsequential limits that a solution to (11.331) exists. Instead, we show that the sequence of iter¬ 
ates actually contracts in certain low-regularity norms, as long as further smallness conditions 
are imposed on 770 and the time interval. This is the content of Theorem 19.61 

Theorems l9.5l and l9.Gl are then combined in Theorem l9.7l to produce the desired high-regularity 
solutions to (11.331) . Indeed, the low-regularity strong convergence and the high-regularity bounds 
combine to show that the sequence of iterates actually converges in a sufficiently high-regularity 
context to pass to the limit in the sequence (without extracting a subsequence) and produce a 
solution to (11.331) satisfying the bounds (12.31) . The high regularity of our solutions requires that 
we impose the compatibility conditions mentioned in Theorem 12.11 
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The transport problem for q 

We study the problem (18.11) in Section [8l the pair (u, p) are assumed to be given, and q is 
solved for. The key feature of the transport equation in a bounded domain is the behavior 
of the normal component of the transport velocity on the boundary: its vanishing allows for 
the construction of solutions without imposing boundary conditions. This vanishing occurs for 
(18.11) precisely because the pair (u,p) satisfy the equations dtp = u ■ M on E, [uj =0 on E_, 
and U- = 0 on 

Using this condition and the high regularity of (u, p), we can produce a solution to (18.11) using 
the method of characteristics. In a C fc -smoothness context we could then simply apply standard 
a priori estimates (which again depend critically on the vanishing of the normal transport 
velocity) to deduce the regularity estimates for q desired in our iteration scheme. However, due 
to the Sobolev regularity ( u,r /), it is not a priori clear that those estimates can be rigorously 
applied to the solution given by characteristics. To get around this technical obstacle, we show 
that the problem (18.11) can be “lifted” to a corresponding problem on the spatial domain T 2 x R, 
and that the restriction of the lifted solution agrees with the solution given by characteristics. 
The advantage of working with the lifted problem is that it is amenable to solution by the 
Friedrichs mollification method, which employs delicate properties of mollification operators that 
are unavailable in bounded domains. The Friedrichs method allows for the rigorous derivation 
of the desired a priori Sobolev estimates, and we deduce in Theorem 18.81 that solutions to (18.11) 
exist and obey the estimates needed for our iteration scheme. 

The two-phase free boundary Lame problem 

The sub-problem <EU takes p as a given function in addition to various forcing terms, and 
produces a solution pair ( u,p ). Notice that (17.11) is not a linear problem: the A coefficients of 
the differential operators and the non-unit normal J\f are both determined by p. We need this 
nonlinearity for two crucial reasons. First, as mentioned above, the equation dtp = u ■ AT is 
essential in using (u,p) to generate the transport velocity in (18.11) . Second, it allows us to take 
advantage of the regularity gain offered by the surface tension terms in (17.11) . Without it, our 
iteration scheme would fail to control p at the highest level of regularity. 

We produce a solution to (17.11) in Section [7] It is simple to see (17.11) as two coupled linear 
problems: a parabolic problem for u with p given, and a transport problem for p with u given. 
This perspective works well without surface tension, and indeed this is how we proceed in 
Theorem llO.il However, with surface tension the regularity demands of the cr±A*r/± terms are 
higher than what can be recovered from the transport equation for p, given the target regularity 
for u. Our way around this difficulty is to recover a solution to (17.11) as a limit as k —>• 0 of the 
k —approximation problem (15.11) . This is accomplished in Theorems 17.11 and 17.21 The former 
establishes the existence of solutions to (15.11) on k— independent intervals, while the latter passes 
to the limit k —>• 0 to solve (17.11) and refines some of the estimates from the former. 

The ^-approximation problem 

Theorems o and 17.21 are predicated on the analysis in Section [5] where solutions to (15.11) 
are constructed, and Section El where k —independent a priori estimates for solutions to m 
are derived. 

Solutions to (15.11) are produced in Theorem 15.11 by way of the contraction mapping principle 
in an appropriate space. The contracting map is constructed by solving the linear sub-problems 
(15.161) and (15.201) . The problem (15.161) is a heat equation; its well-posedness is standard, but we 
record (without proof) the precise form of the estimates of solutions in terms of k in Theorem 
HU of Section 01 The only subtlety in (15.161) is the appearance of the term S. This is inserted 
as a “corrector function” in order to force the compatibility conditions for (15.11) to match those 
of (11.381) for every k € (0,1). The problem (15.201) is studied in Section [3l where solutions are 
produced via a Galerkin scheme and an iteration argument in Theorems 13.21 and 13.41 

The analysis in Section [6] develops a priori estimates for (15.11) that are strong enough to 
allow us to bound the existence time from below in a k— independent manner and to obtain 
k —independent energy and dissipation estimates. Here again the nonlinear structure of m is 
essential in closing the estimates (see for example the handling of the F 3,a terms in Proposition 
16.31) . It is noteworthy that our main estimate, Theorem l6.51 actually fails to provide an estimate 
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of \\u\\ l ^ h an+i, though this term is still guaranteed to be finite. This is due to the fact that 
this estimate depends on being able to simultaneously estimate \\g\\ 2 L oo H 4 ,N+i/ 2 - The structure of 

(15.11) only permits k —dependent estimates of this term, and so we fail to estimate IMli^uv+i 
independently of k. It is only later, in Theorem 17.21 after passing to the limit k —>• 0 that we 
are able to estimate ||?/||^oo^ 4 jv+i /2 and thereby recover an estimate of \\u \\ 2 L 2 H 4 n+i- 

Smallness condition 

Throughout the paper we impose smallness conditions on g$ and the size of the temporal 
interval on which solutions exist. No smallness condition is imposed on uq or q o- The smallness 
of g o is used with a small-time condition to control various norms of g(-,t), which in turn 
provides control of the coefficients of the differential operators in (11.331) (see Lemma IE. II) , a 
Korn inequality (see Proposition IE.3D . and elliptic regularity estimates (see Proposition IE. 41) . 
The smallness of time is used as above and to guarantee that various absorption arguments 
work throughout the paper. 

2.3. Definitions and terminology. We now mention some of the definitions, bits of notation, 
and conventions that we will use throughout the paper. 

Universal constants and polynomials 

Throughout the paper we will refer to generic constants as “universal” if they depend on N, 
fl±, the various parameters of the problem (e.g. g, fi±, a±), p*, p* and the functions p±, with 
the caveat that if the constant depends on a±, then it remains bounded above as either cr± 
tend to 0. For example this allows constants of the form gp + + 3c 2 . + but forbids constants 
of the form 3 + 1 /<t_. Whenever behavior of the latter type appears we will keep track of 
the dependence on a± in our estimates. We make this choice in order to handle the vanishing 
surface tension limit. 

We will employ the notation a < b to mean that a < Cb for a universal constant C > 0. We 
also employ the convention of saying that a polynomial P is a “universal positive polynomial” 
if P(x) = f° r some m > 0, where each Cj > 0 is a universal constant. Universal 

constants and polynomials are allowed to change from one inequality to the next as needed. 

Norms 

We write H k (fl±) with k > 0 and and H S (E±) with s € R for the usual Sobolev spaces. We 
will typically write H° = L 2 . If we write / 6 H k (Jl), the understanding is that / represents the 
pair f± defined on £l± respectively, and that f± € H k (Q±). We employ the same convention 
on Throughout the paper we will refer to the space ofL 1 (fl) defined as follows: 

ofL 1 (fl) = {t)£ fL 1 (fl) | [u] = 0 on E_ and V- = 0 on Eft}. (2-6) 

To avoid notational clutter, we will avoid writing H k (£l) or H k (T ,) in our norms and typically 
write only ||-|| fc , which we actually use to refer to sums 

ll/llfc = ll/+ll/U(Q + ) + II/- ll/U(Q_) or \\f\\k = ll/+ll/U(E + ) + ll/-ll_f/ fe (E_) • (2-7) 

Since we will do this for functions defined on both Q and £, this presents some ambiguity. We 
avoid this by adopting two conventions. First, we assume that functions have natural spaces on 
which they “live.” For example, the functions u, p, q, and fj live on fi, while g lives on £. As 
we proceed in our analysis, we will introduce various auxiliary functions; the spaces they live on 
will always be clear from the context. Second, whenever the norm of a function is computed on 
a space different from the one in which it lives, we will explicitly write the space. This typically 
arises when computing norms of traces onto S-(- of functions that live on fh We use 11 • 11 x to 
denote the norm of the space L p ([0,T]; X). 

Occasionally we will need to refer to the product of a norm of g and a constant that depends 
on ±. To denote this we will write 

7 INI* = 7+ ||»7+||Hfc(s + ) + 7- lk-ll^(s_) ■ ( 2 - 8 ) 

Derivatives 
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We write N = {0,1,2,...} for the collection of non-negative integers. When using space-time 
differential multi-indices, we will write N 1+m = {a = (ao,ai ,... ,a m )} to emphasize that the 
0—index term is related to temporal derivatives. For just spatial derivatives we write N m . For 
a G N 1+m we write d a = d^d^ 1 ■ ■ ■ <9" m . We define the parabolic counting of such multi-indices 
by writing |a| = 2«o + «i + • • • + a m . We will write V*/ for the horizontal gradient of /, i.e. 
V*/ = d\fe\ + <92/e2, while V/ will denote the usual full gradient. 


3. The two-phase Lame problem 

In this section we are concerned with solving the two-phase Lame problem for the velocity 
field u: 


(3.1) 


Here we assume that r/ is given and determines A and AT via (11.211) and (11.261) . respectively, 
and that p is given as well. We also assume that p obeys the estimate 


pdtu — div^4 § = F 1 

in H 

-§ A uAf = F\ 

on 

-ls A u}AT = -Fl 

on 

M = 0 

on 

u- = 0 

on Sb 

u(-,0) = u 0 

in H. 


1 3 

-p* < p(x,t) < -p* for all x € 0 and t G [0, T]. 


(3.2) 


In our analysis of (13.11) we will need to make various assumption about energies associated to 
r) and p. We record these now, first defining a low regularity term that will be used in controlling 
the coefficients of the equations in (13.11) : 

12 


£[V,P\( T ) = SU P 

0 <t<T 


1 1,2 + l|9«l(<)ll?/2 + ||9?»(t)|IL + l|p(i)lll + ll<M<)ll 1 

We also define some high regularity terms: 


(3.3) 


2N+1 

&[v,p](T)= I ( \\v(.t)\\l N+l /2 + I 

j =1 

2 N 


4N-2j+3/2 


+ sup \\r)(t)\\l N + Y\\ d t r l( t ) 

0 <t<T \ 11 


AN-2j+l/2 


and 


£[V,P\(T)= sup 
0 <t<T 



2 

4N—2j 


dt 

2N 

+ sup YW&tPi*) 
2N 


2 

4N—2j 


+ sup \\p{t)\\l N + Y\\ d tP( t ) 
' 1 0<t<T \ 11 


4N-2j+l 


(3.4) 


(3.5) 


We will also need to make certain assumptions on the forcing terms. To describe these we 
first define 

2IV—1 „ 2N 


m) ■■= E I 

3 =0 


4iV—2j— 1 


j =0 


4N—2j—l/2 


and 


2N—1 r 


Foo(t) := E 

3 =0 


d{F\t) 


4N—2j—2 


diF\t ) 


4N—2j—3/2 


We will assume throughout this section that the forcing terms obey the estimate 

rT 


sup Foo(t) + / F 2 (t)dt < 00 . 
0 <t<T Jo 


(3.6) 


(3.7) 


(3.8) 
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We define <9tit(-,0) via 

(pd t u)\ t =o = {dw A §> A u + J F 1 )| t=0 , (3.9) 

and then for j = 0,..., 2 N — 1 we inductively define 

3 - 1 

(pdi +1 u)\ t=0 = - (di~ k pdt +1 u)\ t=0 + dj(div 4 S.AU + F 1 )| i= o. (3.10) 

k =0 


We must assume that ito, i ?1 (-, 0), and F 2 (-, 0) satisfy the following 2 N systems of compatibility 
conditions: 


& t {§AuN)\t= 

=o = &> t F\ |t=o 

on E + 

d{ (S^uAf) 

|t=o = ~diF 2 _ | t=0 

on E_ 

%u 

1=0 = 0 

on E_ 

U- 

1=0 = 0 


on Efc 


(3.11) 


for j = 0,... , 2 N — 1. 


3.1. Weak solutions. Weak solutions will only arise in our analysis in the context of the 
highest-order time derivatives of solutions to m, where they will arise as a byproduct of the 
regularity assumptions on the forcing terms. As such, we will only discuss the meaning of weak 
solution and provide a simple uniqueness result, neglecting a construction of weak solutions. 
We define the time-dependent inner-product on o H 1 according to 

((«>«)) = J J '(*) (| ^A(t) u : B \t) v + l 1 ' div A(t) u div A(t) v) , (3-12) 

where J(t) > 0. Notice that for the sake of notational brevity we have neglected to include t 
in ((it, v)) even though the inner-product actually changes in time. With this definition in hand 
we can give the meaning of weak solutions. 

We say it is a weak solution to (13.11) if 

ue L oo ([0,r];^ 0 )nL 2 ([0,r]; 0 ii' 1 ) and pJd t u 6 L 2 ([0, T]; (o# 1 )*), (3.13) 

u(-, 0) = no, and 

(pJdtu, v)* + ((«,«)) = < JF 1 ,v)* - ( F 2 ,v )_ 1/2 jS (3.14) 

for almost every t e [0, T] and for every v € qH 1 . Here (-, -) m denotes the dual pairing for o H 1 , 
while (•,•)_! ; 2S denotes the dual-pairing on i4 1//2 (S). 

Next we show that weak solutions are unique. 


Proposition 3.1. Weak solutions to (EUD are unique. 


Proof. If u is a weak solution, then we choose v = u as a test function in (13.141) . Employing 
Proposition IE.31 we may easily derive the differential inequality 


4/ n V^ + CM|?</ o 




for some universal C > 0. Then Gronwall implies that 

|2 [ r 


it 


f n P J1 ~Y^ - e ' 


at 


in 


p(0)J(0) 


u -f + J‘\wr. + \p\t /2 


(3.15) 


(3.16) 


for a = C + \\dt(pJ)/(pJ)\\ L oo, with C universal. 

Then if m,it 2 are two weak solutions, then u = u\ — 112 is a weak solution with n(0) = 0, 
F 1 = 0, and F 3 = 0. The equality u\ = 112 then follows from (|3.16l) and the fact that 
min pj >0. □ 
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3.2. Strong solutions. Now we turn to the construction of strong solutions to (13.11) . We say 
that u is a strong solution to (13.11) on [0, T) if 


u G F°°([0, T]-,H 2 n 0 H l ) n F 2 ([0, T]; F 3 ) 
dtu G L°°([0, T];H°) n L 2 ([0, T]-, 0 H l ), 


(3.17) 


u(-,0) = uq, and u satisfies (13.11) almost everywhere (with respect to Lebesgue measure on 
and Hausdorff measure on S-t and £&). 

Our next result establishes the existence of strong solutions. 

Theorem 3.2. Suppose that p satisfies (E2D and that £[rj\(T) < 5 q, where 5o is given by 
Proposition E.f . Suppose that uq G H 2 , F 1 G Lfi’Ek 0 0 L^H 1 , dtF 1 G L^^qH 1 )*, F 2 G 
L™H 1 ! 2 , dtF 2 G L^,H ~ 1 / 2 , and thatu^ andF 2 (-, 0) satisfy the system of compatibility 
conditions 


(3.18) 


Then there exists a unique u that is a strong solution to (EUD- Additionally, dtu is a weak 
solution to 

'pdfidtu) - div.4§4(<9 t t!) = F 1 ’ 1 
-$ A (dtu)M = F 2 ’ 1 
-l§ A (d t u)}M = -Fl’ 1 
{d t uj = 0 
dtu _ = 0 

d t u(-, 0) = [p _1 (div^S^u + F" 1 )] | t=0 in Q, 


— S A uJ\f |f=o — F 2 |t=o 

on E + 

-lS A ujAf \ t=0 = -Fl \ t=0 

on E_ 

N U=o = o 

on E_ 

,u-|t=0 = o 

on E&. 


in 

n 

on 

s+ 

on 

E_ 

on 

E_ 

on 


in 

n, 


(3.19) 


where 

F 1 ’ 1 = dtF 1 - d t pd t u + divg t _4 S A u + divyi S dtA u 
F_|_’ = dtF 2 + E>Q tA uAf + S A udtJ\f 
f!’ 1 = dtF* - [S dtA uj Af - M dtM. 

The solution obeys the estimate 

|2 , ii a n2 , |U t o 2 ll2 


(3.20) 


Ml L°?fr 2 + + 


L™H° + \\dtu\\ L 2 T H 1 + \\pJd t u\\ L ^ oHl ^ 

< (1 + P(K[r,, p](T))) exp(T(l + P(JC[rj, p](T)))) ' 

2112 


2 + Hf 1 " 2 


L™H° 


+ + \\ F2 \\l 2 t H3/2 + \\ dtF2 \\L 2 T H-i/*) » ( 3 ' 21 ) 

where P is a positive universal polynomial such that P( 0) = 0. 


+ \\ F %H^ + \\^ Fl \\l T 


2 (o H 1 )* 


Proof. The proof is a fairly standard application of the Galerkin method. For the sake of brevity 
we will provide only a terse sketch. For more details we refer to the incompressible case (where 
the analysis is actually somewhat more complicated): see for example [11] . |3T|. or [25 j . 

Step 1 - The Galerkin approximation 

Let {wk\kL\ be a basis of qH 1 and for m > 1 set u m (x,t ) = a™(t)w k (x), where the k 
summation runs over k = 1,..., m. Using the theory of linear ODEs we can choose the a7f such 
that 

( pJdtu m ,w k ) 0 + ((u m ,w k )) = (JF 1 , - (F 2 ,w k ) 0 E (3.22) 

and ((u m (0), w k )) = ((uo,w k )) for k = 1 ,m. Here (•, -) 0 denotes the L 2 inner product on 
while (-,-)os denotes the L 2 inner product on £. We then take linear combinations of (13.221) 
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so that Wk may be replaced by u m . This results in the basic energy equality 


7 /* I 77112 /* I 77712 

J / p j2 t~ + + Ur\O 0 - (F 2 ,u"*) os ,. (3.23) 

•7 5T2 l/ n 

Notice that the assumption on £[ 77 ] (T) guarantees that that the smallness conditions of Lemma 
IE. II and Proposition IE.31 are satisfied. Korn’s inequality, Proposition IE.31 allows us to bound 


< 


(( u m ,u m )), 


while Lemma IE. II and (13.211 provide the estimates 


u 


< 


PJ- 


< 


Jn 2 

and trace theory provides the estimate 

-(F 2 , 


u m \\ 2 0 and (JF\u m ) 


u 


) < IlF 2 

/ 0 ,E ~ II 


\U 


< 

0 ~ 


11 • 


,/< 


pj 


\u 


(3.24) 


(3.25) 


(3.26) 


We combine (I3.23|1(I3.26|) with Cauchy’s inequality to absorb the ||'u m j | 1 term onto the left; 
applying Gronwall’s lemma to the resulting inequality results in the bound 


u 


|2 
I L° 


>ho + h m \\ 2 L * T Hi Z exp(T(l + P(K[ V ,p](T)))) (llrxollo + ll^ll^o + 


\l 2 t h° 


(3.27) 

Next we differentiate (13.2211 in time and use a similar argument to get estimates of the sum 
\\dtu m \\ L oo H0 + \\dtu m \\ L 2 H i. Two complications arise. The first is that differentiating (13.221) 
introduces a number of commutators. These may be dealt with via the estimates (13.2711 and 
the usual Sobolev embeddings. The second is that we must estimate ||<9tTi m (0)||Q in terms of uq. 
This can be accomplished by using (|3.22ll . evaluated at t = 0, in conjunction with the choice 
of initial condition for u m { 0) and the compatibility conditions (|3.18l) . The resulting estimate 
bounds \\dtu m \\ 2 L oo H o + || d t u m " 2 


\l 2 h 1 by the right side of (13.2111 . 


Step 2 - Passing to the limit and higher regularity. 

The estimates from Step 1 are uniform in m and hence allow for the extraction of weak and 

o-i 


weak-* limits. Passing to the limit yields a function such that u,dtu € L^?iL u n L 2 nH l and 


(pJd t u, w ) 0 + (( u,w )) = (JF 1 ,w) q - ( F 2 ,w ) 0 J 


(3.28) 


for a.e. t € [0,T] and every w € qH 1 . The L^H 0 and LrpH 1 norms of u and dtu are controlled 
by the right side of (13.2111 by virtue of lower semicontinuity. We then apply the elliptic regularity 
result for the Lame problem (a variant of Proposition IE.41 adapted to the boundary conditions 
involving §, which can be proved in a similar manner by appealing to Section 3 of m) to deduce 
that 11 ^ 11 ^ 00^2 + \\u \\ 2 L 2 H 3 is also bounded by the right side of (13.211) . The only term remaining 

to estimate in (13.211) is ||pJ9 2 u ||^ 2 This can be accomplished by using ([3.2811 and the 

existing estimates. The equation (I3.19H is seen to hold by differentiating (13.2811 in time. □ 


3.3. Higher regularity. Our aim now is to show that the strong solutions of Theorem 13.21 are 
of higher regularity when the forcing terms are more regular. We first define the forcing terms 


that appear for the time differentiated versions of (13.11) . We set 

F 1 ’ 0 = F 1 , F 2 ’ 0 = F 2 , and F 2 ’ 0 = F 2 , (3.29) 

and then for j = 1,, 2 N we iteratively define 

F lj = 5 t F lj_1 - d t pd{u + divQtA^Adi^u + div^(3.30) 
Fl’ j = d t F 2,j ~ 1 + S dtA di~ l uM + S A d{- l ud t N (3.31) 

F 2 J = dtFl*- 1 - Ifsa^ar 1 J M - Ik^rJ dtM. (3.32) 


The following lemma will be used in conjunction with an iteration argument to improve the 
regularity of strong solutions to m when the forcing terms are of high regularity. It is a 
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simple variant of Lemma 3.2 of [25], which is a technical refinement of Lemma 4.5 of [ 11 ] . As 
such, we omit the proof. 

Lemma 3.3. For m = 1, ..., 2N — 1 and j = 1, ..., m, we have the following estimates, where 
P is a positive universal polynomial such that P( 0) = 0: 


|i^|| 


L 2 iH 2m-2j+l + \\F J |lz,2 jj41V-2j+3/2 <(1 + P(%,p](T))) 


( r T 


2 ^ || 

2 


2 \ 

[J P 2 {t)dt + 

d{u 

+y: dfu 

L 2 ^2m-2i + l || 1 

T £=0 

L|,_ff2m-2j+3 + 


L?yH 2 rn ~ 2 i + 2 J 


(3.33) 


I pl,j\\ 2 II rp2,j II ^ 

I II L™H 2m ~ 2 j “HI Wl^h 2 ™- 2 ^ 1 / 2 


(1 + P(A[r,,p](T))) sup -Foo(f) + 

\ 0 <t<T 


d^u 


j -1 


L oo^2m-2j 


+E p 


u 


£=0 


L™H 2rn ~ 2 i+ 2 


(3.34) 


and if we denote the functional v eA (JF 1 ,m ,u ) 0 — (F 3 ,m ,u) s in (o H 1 )* by F m , then 


d t F m 


< 




(i+p(%, P ](T))) ( [ T p 2 (t)dt+Y^ h 

V° £=0 11 


u 


LJ?H 2 


+ 


din 


L^H 3 


+ WU\\ 2 L 2 TH 2 + \W +1 U" 2 


L^HO 


(3.35) 


We can now state our higher regularity result. 

Theorem 3.4. Assume that uq £E H iN , that F 1 ,F 2 satisfy (13.81) . and that the compatibility 
conditions (13.111) are satisfied. Assume that p satisfies (13.21) . and that £[ 77 ](T) < do, where do is 
given by Proposition \ E.fi Then there exists a unique u satisfying dfu € L°°([0, T]; H iN ~^) n 
L 2 ([0,T]-,H 4 N ~ 2 i +1 ) for j = 0,...,2JV and pJd 2N+1 u € L 2 ([0, T]; ( 0 H 1 )*) so that dj/u solves 

pdt(d{u ) — div^S^c^n) = F l)J in Q 
-§ A (d^u)Af = F 2J on E+ 


§A{dt u ) 
= 0 


AT = -F 2 ’ j on E_ 

on Li¬ 
on St 

strongly for j = 0,..., 2 N — 1 and weakly for j = 2N. The solution obeys the estimates 

2 

L 2 T ( 0 my 


(3.36) 


2 TV 


E K 

3 =0 


u 


LfH™- 2J+1 


+ 


dju 


dlu _ = 0 


pJdj N+1 u 


2N 


< 

r^j 


(l + P(%,p](T)))exp((l + P(g[ 7 ? , / o](T)))T) ( £||9t'<,0) 

i=o 


+ 

4N-2j Jo 


J P 2 {t)dt 


(3.37) 


and 
2 TV 

\\ d t u 

j =0 


L oo H4 JV-2j 


< (1 + P(J% p]{T))) exp ((1 + P(<£[ 77 , p](T)))T) 


2N 

,1=0 

where P is a universal positive polynomial such that P( 0) = 0 


+ sup Poc(f) + [ P 2 (t)dt I , (3.38) 

4A r -2j 0 <t<T Jo 
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Proof. The proof follows the same basic outline (though the arguments here are somewhat 
simpler), as that of Theorem 4.8 of |11| : we iteratively apply Theorem 13.21 Lemma 13.31 to 
estimate the forcing terms appearing for each j, and the elliptic regularity result. We omit 
further details for the sake of brevity. □ 


4. Heat estimates 


Suppose that 0 < k < 1. In this section we are concerned with the problem 


d t p — kA*77 
= Vo 


where po G H 4JV+1 (S) and / satisfies 


/ in £ x (0, T) 
in £ 



2 

dt < oo. 

4N—2j 


(4.1) 


(4.2) 


A simple interpolation argument (see for example Lemma A.4 of El) allows us to deduce from 
(14.2|) that dif G C°([0, T]; H AN ~ 2 i~ l {YP)) for j = 0,..., 2IV — 1. We may use this in (14.111 to 
inductively define dfp(-, 0) G H 4JV_2 - ?+1 (£) for j = 1..... 2N 

We now state a result on the well-posedness of (14.11) and estimates of its solutions. The proof 
is standard and thus omitted. 


Theorem 4.1. Suppose that po G H 4N+1 (T,), that f satisfies (14.211 . and that the data dtp(-,0) G 
//4V-2.7+1 (£) are determined as above for j = l,...,2iV. The problem (14.111 admits a unique 
solution that achieves the initial data &[p(-, 0) for j = 0,..., 2 N and that satisfies the estimate 



2 

L oo H 4JV-2i 


+ 


2IV+1 

£ 

3 = 1 


dj/n 


2 

1,2 J/4JV-2J+2 



' 2 N 

2 2N 

2 

+ K, 

Y ^ 

3= 0 

+ Y W&iv 

LSf H 4N ~2j+i ^ || * 

T 3=0 

L 2 tH 4N-2j + 1 


+ WvW 2 L? r H iN + 2 


< 


2 TV 


3=0 



(4.3) 


5. The k —approximation problem 


Our goal in this section is to produce a solution pair (u, p) to the problem 


' pdtu — div _4 u = F 1 

in H 

dtp — k,A*t/ = u ■ AT — kE 

on £ 

—S_AuJ\f = —a + A*pJ\f + F 2 

on £ + 

— A/ - = <j—A*pAT — Ff 

on £_ 

N =o 

on £_ 

u- = 0 

on £fe 

M',0) = u 0 ,p(-, 0) = p 0 . 



(5.1) 


Here we assume that uq G H lN {Pl), po G H in+1 (Y<), F 1 , F 2 satisfy (13. 811 . Notice that the 
problem (15.11) is nonlinear: the terms A and JV are viewed as being generated by p itself. We 
will again assume that p is a given function on the given time interval [0, T*] and satisfies (13.211 
(with T replaced by T*) in addition to the bound 


<£[p](T*) := sup £oo[p(t)] 
o <t<T* 


sup 

0 <t<T, 




4N-2j+l 


<l + P(Eg). (5.2) 
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We now show how to construct the data 0), c 0)) for j = 1,, 2N from the pair 

(rt 0 ,r? 0 ) in a manner that is consistent with the compatibility conditions for the problem (11.331) . 
Given (<^?rt(-, 0), 0)) for some j = 0,..., 2N — 1 we first use the first equation in (15.11) to 

solve for c^ +1 u(-,0). Then we set 

dt +1 v(-, 0) =d}(u- AT)\ t =o- (5.3) 


In this way we iteratively define all the time-differentiated data for j = 1,... ,2 N. However, 
the only way this can be consistent with the natural compatibility conditions needed to study 
(15.11) is if 5 : £ x [0, oo) is chosen so that 


S(.,0) 

8iS(; 0 ) 


A*r/o G H 4 1 (E), and 


(^dj l (u • AT)) ^ € H 4N ~ 2J-i/2(e) for j = 1,..., 2 N. 


According to Proposition IF.31 there exists a S satisfying (15.41) as well as the bounds 


(5.4) 


21V—1 2 2N r 

Y sup dfs(t) +y: / 

^ t> 0 4N-2j-l ^ In 

j =0 >u j=0 J{J 


dim 


4N-2 j 


dt 


2N 


< ||A*r/ 0 ||4 Ar _ 1 + Y || A * y d t 1 ( u ■ Y)J (0) 

3 = 1 


41V—2j — 1 


(5.5) 


In order to produce a high regularity solution to m we must assume that the uo, r/o, and 
the forcing terms F l , F 2 satisfy the compatibility conditions 


” -di(§> A uJV)\ t=0 = di (—o- + A*r/Af + Ff) | t=0 on £ + 
d{ (S auN) |i=o = d{ (a-A*r]Af - F 2 ) | t=0 on £_ 
U=o = 0 on £_ 

i=o = 0 on 


dt u 


K&tU- 


(5.6) 


for j = 0 ,..., 2N — 1. 
For 0 < T we define 


r T ( 2N ii 

2i 2 m(t)=/ qZ Mm 


G=° 

2N 


2too[u](T) = sup Y\\dim 


0<t<T J =0 


4N-2j+l 
2 

4N—2j 


+ 


pJd? N+1 u(t ) 


(off 1 )* 


dt 


(5.7) 


Similarly, we write 

f.j' 2JV+1 


T 2N+l 2JV 2 

W(T) = / <9 t J ??(f) dt and ©oofoKT) = sup ^ fe(t) ■ (5.8) 

40 1 41V-2.J+2 0<t<T “j) II 4IV-2J+1 


We then define 

2l[u](r) = 2l 2 M(T) + 2tooM(T) and ®[r/](T) = <B 2 M(T) + ®ooW(T). 
Finally, we also define 


21V—1 


21V—1 


°^ T j=0 


5M(T) = sup ^ \ d 3 t u{t) and 3[r]](T) = sup Y W) 

n<t<T I 41V—2j —1 n<KT I 


J=0 


41V-2j-l 


(5.9) 


(5.10) 
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5.1. Constructing the mapping. We will employ a fixed point argument in order to produce a 
solution to the nonlinear k —approximation problem (15.11) . We thus begin with the development 
of an appropriate metric space in which to work. 

Notice that if ^MCO+^CCK? 1 ) < °o, then the maps [0,T] 9 t ha dfv(-, t) € H 4N ~ 2] (Q) and 
[0, T] 9 t i-A dlC(-,t) € H 4 N ~ 2 ^(E) are absolutely continuous for j = 0,... , 2/V. In particular, 
0) and <9^ £(•, 0) are well-defined for j = 0,..., 2 N. This then allows us to define the metric 
space, for 0 < T < T*, M\ , M 2 > 0 and r > 2 an integer, 


X(M 1 ,M 2 ,r,T) = {(v,C) | 21M(T) < M^,®[C](T) < M 2 tW 2 , 

and d}v(-, 0) = ^u(-,0),^C(-,0) = ^r?(-,0) for j = 0,... ,21V}. ( 5 . 11 ) 


The space £(M\, M 2 , r, T) is complete due to the fact that {(u,C) I 21M(T) + 5S[C](T) < 00 } is 
a Banach space. 

We now define a mapping M : X(Mi, M 2 ,r, T ) -A- X(Mi,M 2 ,r, T) as M.(v, () = (it, 77 ), where 
77 and u are determined through the following two steps. In these steps we will assume that 
0 < T < 1 and that £[ 170 ] < 5q/2, where 5q is given by Proposition IE.41 
Step 1 The rj equation 
Note that 


2 TV p 1 2 N p 

E/ a?(w-JV[CD(0 * _ rft<(i + <MC](T))E / SM*) 


J=0" u J 1=0' 

The usual Sobolev interpolation allows us to estimate 


47V-21+1/2 


di. (5.12) 


d(v(t) 

2 

< 



dfv(t) 


4JV-21+1/2 ~ 


AN—2j 



4IV-21+1 


when j = 0,..., 2IV, and hence 


2A „ T 

§1 


^u(t) 


2iV 


dt ^i 


AN-23+1/2 —j o 


2iV -j 1 


1/2 


< 


f 9M*) 

U=0 Jo 


AN—2j 


dt 


AN—2j 

2N 

E 

l=o 


di v (t) 

n -lO 


4JV-2/+1 


dt 


1/2 


dt 


(5.13) 


4JV-21+1 

<VTvCmV@). (5.14) 


Combining the estimates (15.121) and (15.141) . we find that 


2N 


1=O' 


E/ «?(” -vicdw 


A < (i + <B„[<](r))\A ; v'2ioo[t.](r)v'a 2 H(r). 

47V-2j 


(5.15) 


Estimate (15.151) and (15.51) allow us to use Theorem 14.II with f = v • A/"[C] — and 5 defined 
as in (15.41) (15.51) to produce a unique solution 77 to 


dtr] — kA*t] = v ■ 1V"[C] — kE in E x (0, T ) 
??(•, 0) = 770 in E. 


(5.16) 


The theorem guarantees that c^ + 1 ? 7 (-, 0 ) = d{(v ■ AT)|t=o = d{ (it • AT)|t = o for j = 0 ,..., 2 N — 1 , 
with the understanding that this last quantity is the initial data computed previously in (|5.3|) . 
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and that rj obeys the estimate (utilizing (15.151) 1 


2N 

k 2 (® 2 [r/](T) + ( Boo[r?](T)) < e T £ ||«-, 0) 

1=0 

.j 1 2AI 


4AT-2j+l 


1=0 


4IV-2 j 


2AT-1 

4IV+1 + | dt( u ' AO I 

1=0 


< 

r^j 


+ [ e T *^te(u-AT[C]-«S)(t) 

J 0 : n 11 

+ (1 + ® oo [C] (T)) vT VSloo M (T) M (T) 


4AT-2J-1 


2IV-1 

4JV+1 + Y1 | d t( u 'N) | 
1=0 


4IV—2j—1 


+ 1 + 


where in the last line we have used the bounds built into the definition of 3L(M \, M 2 , r, T) and 
the assumption that T < 1. 

From (15.171) we see that if M 2 is sufficiently large (depending on the data, < 7 , etc) and T is 
sufficiently small (depending on k and Mi, M 2 ), then 

®wcn<^( 5 . 18 ) 

ru 

which is the requirement for tj to be part of a pair belonging in 3L{M \, M 2l r, T). On the other 
hand, from Lemma lD.il 


2[V\(T) < £[r? 0 ] + [ T &oo[ri(t)}dt < £[%] + < ^ + ^ < <5 0 

jo 2 

if T is further restricted. 

Step 2 - The u equation 

Now we seek to solve the problem 


pd t u - div^ [r?] S A[v] u = F 1 

in D 

-§> A[v] uJ\f[r]} = -cr + A*riAf[r]\ + F% 

on 

~ NM = cr-A„r]Af[T]} - F 2 

on £_ 

H = 0 

on 

u- = 0 

on Tb 

u(-,0) = u 0 , 



(5.19) 


(5.20) 


where rj is the function constructed in Step 1 and A/" [ 17 ], -4. [ 17 ] are determined in terms of rj as in 
(11.211) and (11.261) . Let us write F 1 = F 1 , F 2 = — <j + A*r]J\f[r]} + F 2 , and F ’ 2 = a-A*r]J\f[r]\ — F 2 . 
Then we write F 2 (t) and (t) as in (13.61) and (13.71) . except with F 1 and F 2 in place of F 1 , F 2 . 
Then clearly 


2N 


2 N 


A'oo (t) < Foo(t) +1 + 5^ mv(t) 


1=0 


4IV-2j+l 


Mvit) 


1=0 


4N-2j+l 


2N 


£2 it) < F 2 (t) + (1 + C 7 2 ) ( 1 + || %v(t) 

1=0 


4AT-2J+1 


2 N 

1=0 


(5.21) 


4IV—2J+3/2 
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A simple interpolation argument, in conjunction with the Cauchy-Schwarz inequality, pro¬ 
vides us with the estimate 


r T 2N 


3=0 


4N-2j+3/2 


r T 2N 

dt ~ I 


/ 0 


3=0 


4N-2j+l 


dlv(t ) 


4N-2j+2 


dt 


< < Vtvb 2 [ v \(t) + ©ooW(r)). ( 5 . 22 ) 

We then deduce from (15.211) and (15.221) that 

(5.23) 


sup (t) < sup Too(t) + (1 + 03 oo [?7](r))Q3 oo [77](T) 

0 <t<T 0 <t<T 


and 


[ Mt)dt < [ T 2 (t)dt + (1 + a 2 )(l + ® 00 [ ?7 ](r))\/f(®2[7/](r) + ® 00 [j 7 ](r)). (5.24) 
Jo Jo 

Now, owing to the bounds (15.231) - (15.241) . the compatibility conditions (15.61) . and the estimate 
(15.191) . we may apply Theorem 13.41 with F 1 and F 2 in place of F l ,F 2 . This yields a unique 
solution u to (15.201) satisfying (after combining (13.371) (13.381) with (15.231) (15.241) 1 


Mu}(T) < (1 + P(®[r?](T) + <Z\p](T))) exp (T(l + P^MT) + £[p](T)))) 

2N 9 rT 

J2 N“( - > 0 ) 


X 


J=0 


2 

4N—2j 


+ J F 2 {t)dt + (1 + a 2 )VT<B[r]\(T) J , (5.25) 


and 


SW«]CO < (1 + Pmv\(T) + <£[p](T))) exp (T(l + P&MiT) + (£[p](T)))) 


x 


2 N 

J=0 


+ sup P 0 O (t) + m 00 [ri}(T)+ / P 2 (t)dt + (1 + a 2 )Vr<3[r,\(T) 
4 N—2j o < t <T Jo 


(5.26) 

for some universal positive polynomial P(-) with P(0) = 0. 

Recall that we assume the bound Then by restricting T further and employing the 

estimate (15.181) . we may deduce from (15.251) and (15.261) that there exists a natural number i 
such that 


2l[u](T) 


1 + P(£S) + K" 


2 N 


<3=0 


2 

4N—2j 


+ sup Foo(t) + [ p 2 (t)dt + 1 j . 
i 0 <t<T Jo 


If M\ is sufficiently large and r is set to £, then 


Mi 

*[«](T) < -T: 

ru 


(5.27) 


(5.28) 


which means that u satisfies the estimate required to belong to a pair in X(M\. M 2 . r. T). 

Step 3 - Inclusion in 3£(Mi,M 2 ,r,T). 

In particular, (15.181) and (15.281) imply that if Mi, M 2 > 0 are taken to be sufficiently large, and 
0 < T < 1 is taken to be sufficiently small, then (it, 77 ) G X(Mi, M 2 , r, T). This allows us to set 
M(v, C) = (u, rf) G X(Mi, M 2 ,r, T), and so the mapping A4 : £(Mi, M 2 ,r. T) —y X(Mi, M 2 ,r, T ) 
is well-defined. 

5.2. Fixed point. Now we show that the mapping constructed above has a fixed point. 

Theorem 5.1. Assume that p is given on the time interval [0, T*] and satisfies (13.21) and (15.21) . 
Let <5o be given by Proposition E.fi and assume that 6 G (0, do). Suppose that that £[ 770 ] < 5/2, 
where £ is given by (IA.10D . There exist M\,M 2 > 0, r G N, and To = Tq{k,g,5,£q) G (0,T*] 
such that if 0 < T < Tq then M : X(M\, M 2 ,r,T) —> 3L(M\, M 2 ,r,T) is a contraction, and 
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therefore admits a unique fixed point ( u,rj ) = M{u,p). In particular, there exists a unique pair 
(u, rf) satisfying 

21 [u](T) < M\K~ r and [ 77 ](T) < M 2 k ~ 2 (5.29) 

that solve (EH and achieve the initial data cP t u(-, 0) and cP t ri(-, 0) for j = 0,, 2 N. Moreover, 

£[v]( T ) < $ and 3[j?](T) + 3MCO < £°[u 0 ,vo] + 1, (5.30) 

where £°[uo,po] is as defined in (IA. 131) . 


Proof. Let (uj,Ci) G £(Mi, M 2 , r, T), ( Ui,rji ) = M{vi,Ci) G X(M 1 ,M 2 ,r,T) for i = 1,2. Then 
u := u\ — u 2 and 77 := 771 — 772 satisfy 


and 


pd t u - div^[ m] = il 1 in 0 

= -a + A*r/AA[r/i] + on E + 

- [S^p] AAfri] = (T„A*r/AA[? 7 i] + il 2 on E_ 

[it] =0 on 

U- = 0 on E& 

pM) = o 


dtp - kA^p = (vi - u 2 ) • Af[Ci] + v 2 ■ (AA[Ci] - AA[C 2 ]) 

7?(') 0) = 0, 


(5.31) 


(5.32) 


where 

H l = div A[m] _ A[m] S^[r, 2 ]U 2 + dW A[m] §A[vi}-A[r 1 2 } U 2 

H'l = -cr + A*r/ 2 (A[7?i] - M[p 2 \) - S^ [?7i] 77 2 (AA[t7i] - AA[r/ 2 ]) - S A[m] _ A[m] u 2 J\f[p 2 ] (5.33) 

#- = o'-A*7 ?2 (A[ryi] - AA[r? 2 ]) - [S^ [rjl] « 2 ] • (J\f[p 1 ] - A f[p 2 }) - [S^]-.^]^] Af[p 2 }. 


Notice that the H terms cancel in (15.3211 since they are both computed from the same initial 
data. 

From (15.151) in Step 1 above, we know that 


2N rT 

E / ^((«i-«)-JV[Ci])(<) 


j= 0 J0 


47V —27 


dt 


< (1 + ® oo [Cl](F))\/T 1 ^/2loo[l7i - u 2 ](r)V2l 2 [wi - V2\(T) 

< (l + v / T2l[r; 1 - u 2 ](T). (5.34) 


A similar argument shows that 
27V „ T 

E/ ^(wMCi])-y[c 2 ])(() 


47V—2j 


dt 


< ©ooICi - C 2 ](T)x/fVSlooM(T)v/2l 2 [u 2 ](T) < -AVTSooICr - C 2 ](T). (5.35) 

Hi 

Then, as in (15.171) . we have the estimate 

®2fo](T) + ®«MT) < - v 2 ](T) + - c 2 ](r). (5.36) 

\Hi Hi J Hi 

From (|5.25p in Step 2 above, we know that 


2l 2 M(T) < (1 + P(®fo](T) + (£[p](T))) exp (T(l + P(*B 00 H(T) + £[p](T)))) 

x(jf H 2 (t)dt + (l + a 2 )v/T s B[r ? ](r)^ (5.37) 
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for some universal positive polynomial P with P(0) = 0, where we have written 
2JV—1 0 2N 


2 

4N—2j — 


n 2 (t)= E I %H\t) 

3=0 

Similarly, (15.2511 in Step 2 yields 


i + H^ i (t)|| (o ^+E||^ 2 w 


3=0 


4N—2j—l/2 


2tooM(T) < (1 + P(2%](T) + £[p](T))) exp (T(l + P(®oo[r?](T) + £[p](T)))) 


where we have written 


2N-1 r 


P-oo(t) E! 

3=0 




4N—2j—2 


+ 


djH\t) 


4N—2j—3/2 


(5.38) 


x f sup Hoo(t)+ H 2 (t)dt + , Boo[i?](r) + (l + ( T 2 )v/T53[r ? ](r) ) , (5.39) 

yo<t<T Jo J 


(5.40) 


We now seek to estimate the % terms appearing on the right side of (15.3711 and (15.3911 . 
Standard nonlinear estimates lead us to the bounds 


n 2 (t)dt < ®ooM(T)(i + P^oobnKT) + ®ooM(t))) 21 2 [u 2 ](t) + ®ooM(t)® 2 M(t) 


< p 


M 2 M 1 M2 Mi 

K r k'^-I- 2 i..2 


®ooW(r) (5.4i) 


and 


SUP H^t) < ®ooW(r)(l + ®oohl](r) + ®oo [l? 2 ] (P))2loo [u 2 ] (T) + ^ [ V ] (T)®*, [t? 2 ] (T) 
0<t<T 


< P 

r-^j 


(M 2 M\M 2 M\ \ 

+ K r + 2 + "k- 2 / 


®ooN(T). 


(5.42) 


Now we sum the estimates (15.371) and (15.391) and then combine the resulting estimate with 
(15.41 p and (|5.42D to deduce that 


21 M(T) 


, M 2 
1 + P — £ 


P 


exp 

M 2 


T (1 + P(^ 


M1M2 Mi 

K r Av r + 2 K 2 


$ B 00 [i?](T) + (l + u 2 )%/f®[? 7 ](r) 


(5.43) 


for some universal positive polynomial with P(0) = 0. Combining (15.3611 and (15.431) then yields 
the estimate 


2tM(T) + ®W(T) < (l + P (^)) exp (r (l + P (^|))) 


x 

x 


(2f + %■) ~ «]cn + 




(5.44) 


Finally, from (|5.44|) we see that if we further restrict T in terms of Mi, M 2 , l, &, and k, then 
2l[ui - u 2 }(T) + <B[ m - 772 ](T) < \ (2t[ui - v 2 ](T) + ®[Ci - c 2 ]CO), (5.45) 

which implies that the mapping M : X(Mi, M 2 , r, T) — > 3 L{M \, M 2 , r, T) is a contraction. 

The existence of a fixed point satisfying (I5.29|) is then an easy consequence. The estimates 
(I5.30P follow from (I5.29p . Lemma lD.il and standard estimates of the data. □ 
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6. Estimates for the k —problem ()5.1f) 

Our goal in this section is to derive k —independent estimates for the problem (15.11) . These 
will eventually allow us to pass to the limit k —>• 0. In this section and the next we must work 
with a slightly weaker form of the dissipation for it, which is defined as 


V[u\ = \\u\\ 4N + 




2 N 

3 = 1 


4N-2j+l 


Here and at several points in this section we employ the notational convention 

12 iiqa„.i|2 


= X) ll 8 ”“l 


aeN 2 

\a\<m 


( 6 . 1 ) 


( 6 . 2 ) 


for any m > 0 and any norm ||-||. 

Recall that 5 is the function defined in (15.41) (15.51) in terms of the data. We will need to refer 
to the following functional: 


2N-1 .2 2N fQO 

3 ■■= sup 3s S (*) +S / 3s 3 (*) 

U t>0 47V_2j ” 1 U Jo 


AN-2j 


dt. 


It’s easy to see from (15.51) that 


3< ^\\V°V*rio\\l N + P{£S[u 0 ,rio]) 


for a positive universal polynomial P such that P(0) = 0. 


(6.3) 


(6.4) 


6.1. Preliminaries. Rather than work directly with the solutions from Theorem 15.11 we will 
prove our estimates in a somewhat more general context. We assume the following for some 


parameter 5 > 0 and time interval [0, T*]. 

The parameter 6 > 0 satisfies 5 < 5q where Jo is given by Proposition IP. 41 (6-5) 

The initial data 0) and 0) for j = (),..., 2N are given as before and ( 6 - 6 ) 

satisfy the compatibility conditions (15.61) . 

The data satisfy £[%] < <5/2, where £[ 770 ] is given by (1A.10D . (6-7) 

A solution (u, 77 ) to (15.11) exists on the interval [0, T*] and achieves the initial data. ( 6 . 8 ) 
The solution satisfies the estimate 2l[u](T*) + ®[t 7 ](T*) < 00 , where 21 [it] and 25[ 77 ] 

are as defined in ([5.91) . (6-9) 

We have the estimates £[ 77 ](T*) < 5 and 3 [? 7 ](T*) +3[ti](T*) < P(£ q) + 1 
for a universal positive polynomial P such that P(0) = 0. Here 3[ 77 ],3[it] 
are defined by (15.101) . (6.10) 

r T * 

The forcing terms satisfy sup Poo(t) + / p 2 (t)dt < 1 + P{£q) 

0 <t<T* Jo 

for a universal positive polynomial P such that P(0) = 0. (6-H) 

The function p is defined on [0, T*] and satisfies m an d m- (6-12) 


The assumption (16.91) guarantees that (it, 77 ) are regular enough for us to apply d a to (15.11) 
for a € N 1+2 with |a| < 4 N. This results in 
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pdtd a u — div 4 §_4 d a u = d a F 1 + F 1,a in Q 

dtd a r] — kA= ( d a u ) • M — nd a E + F 3,a on E 
-S A d a uAf = —cr + A*d a r]Af + d a Fl + F+ a on E+ 
- lS A d a uj A f = a-A*d a r]Ar - d a Ff - F 2 ’ a on E_ 
[<9 Q u]] =0 on E_ 

d a U- =0 on Eft, 

where for i = 1, 2,3, 

F- ^ C a fid a ~Ppd t d p u,i + ^ C a ,p9 a ~ 3 A( m d 13 d m ( pA( k df.Ui ) 


(6.13) 


/3<a 

+ (// + ^ + d a ^dk(A£ m d m U£)^ , (6.14) 

/3<a 


P<a 


F i“ = Y C a , p ( M 9 a -^(A^Afe)^4^ + iJ^-^MtAi^dPdkUi 


/3<a 


+ E C «,/3 ( (V - y) F-tWiAtd&dkUt - cr+cP^A/^A*^ (6.15) 


and 


F - a i = E ( d a - p WAk) +5 q -^(A4Ya ; ) L^^fcUi' 


P<a 


+ Y C ( d a - p {M % Aik) 

/3<a 


p' ~ y- ) d^dkiii 


+ (T-& x - p Af i d p A^. (6.16) 


Also, 


F 3 ’ a = Y C a ,pd a ~PNd? 

P<a 


U. 


(6.17) 


6.2. The estimates. We begin with a basic energy identity. 

Lemma 6.1. Assume ()6.5f) (16.121) . Let a € N 1+2 satisfy |a| < 4 N. Then 


d ( f pJ 


dt \Jn 2 

+ K 


I^I 2 + /J |^' 2 + ||V^' 


2 ' + [ \lD A d a u \ 2 + Jp,'\div A d a u \ 2 

J n 2 


2 , _ i a oa„|2 / dt(pj) pJ ,„ a ,2 


[ \V*d a p \ 2 + a\A*d a ri \ 2 = I 

Jt, Jn 


\d a u\ z + Jd a F L -d a u- / d a F z -d a u 

n Pd 2 ,/ s 


+ [ JF 1,a ■ d a u — [ F 2 ’ a -d a u + [ d a pd a u- A7+ [ (nd a E - F 3 ’ a ){-d a p + aA^rf). 
Jn Jt, Jt, J e 

(6.18) 


Proof. We multiply the first equality in ()6.131) by Jd a u and integrate over fi. After integrating 
by parts and using the boundary conditions, we find that 

^ (./ l^ Q “| 2 + J \ + /^|DPW + V|div^«| 2 + i«/a|W 

= [ ^^^-\d a u \ 2 + Jd aFl ■ d a u- [ 8 a F 2 -d a u 

Jn pJ 2 J E 

+ f JF 1,a ■ d a u — [ F 2 ’ a -d a u+ [ (Kd a ~ - F 3 ’ a )crA*d a rj. (6.19) 
Jn J e 4s 
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Also, we multiply the second equality in (16.131) by d a rj and integrate over X. After again 
integrating by parts, we find that 

4- I h da 'n\ 2 + fi [ N*d a v\ 2 = [ d a r 1 (d a u-N) + (-Kd a Z + F 3 ’ a )d a r ] . (6.20) 
dt Jy 2 Jy Jy 

Summing (16.191) and (16.201) yields ()6.181) . □ 

Our next result provides some estimates of the forcing terms that appear in the equations 
(16.131) as a result of commutators with d a . The proof is similar to that of Lemma 13.31 and the 
H forcing estimates of Theorem 15.11 and is similarly omitted. We recall that £oo[p\ is defined 
in (15.21) . 

Lemma 6.2. Let a G N 1+2 satisfy |a| < 4 N and let F l,a be given by (16.141) - (16.171) . Then there 
exists a polynomial P with universal positive coefficients and P( 0) = 0 such that 

IT'llo + lh 2 ’”l| 2 -i/2 S (r(£>]) + H + £>]) + mwb [»]• (6.21) 

Additionally, 

|| (j p3,a + (^v.tj) • u)|| 2 < S°[rj\S[u] + 2[ V }V[u), (6.22) 

and 

||V*(F 3 ’“ + (d Q V*r?) • «)||5 < + £[rj\T>[u]. (6.23) 

With this lemma in hand, we can turn the energy identity of Lemma f6.1l into a useful estimate. 
Note that in this proposition we employ the notation defined in (16.21) . 

Proposition 6.3. Assume (16.51) - (16.121) . Then for 0 < T < T* we have the estimate 


2 N 

3=0 


U 


LfH° 


2N 


+ 




LfH 1 


+ 


di/n 


L oo H 4JV-2j 


+ a 


dj'V^r] 


L oo i?4 AT-2j 


3=0 


k j v*r 7 

2 

+ <T 

<9* A*r? 

11 

L 2 t H 4JV-2J 





. ^ e /T £f[u 0 l Vo] 


+ e^ T / {F 2 (t) + (P(£°[ri)) + €oo \PWH + S 0 [rj\) + P(£[rj\)V[u]) dt, (6.24) 

Jo 


where 7 = C (1 + J[rj\(T) + l£[p](T) + ^[u](T))), and P is a polynomial with positive universal 
coefficients satisfying P( 0) = 0. 


Proof. Let a G N 1+2 satisfy |a| < AN. Taking (16.181) of Lemma l6.ll as our starting point, we 
seek to estimate each term on the right hand side and to estimate non-time-derivative term on 
the left. In our subsequent analysis we will rewrite F 3,a = (T 3 ’" + (< 9 “V* 7 ) • u) — (cPV*? 7 ) • u. 
First note that we may use Proposition IE. 31 to bound 

||3“u||?< [ ^-\B° A d a u\ 2 + Jp'\div A d a u \ 2 . (6.25) 

Jn 2 

This will allow us to absorb most of the d a u terms appearing on the right side of (16.181) . Using 
the definition of P 2 , given in (13.61) . in conjunction with the estimates (16.211) and (|6.22l) of Lemma 
16.21 and the bound 0 < n < 1 , we may bound 


+ \\d a F 2 \\ 2 0 + ||F 1 ’“||g + ||T 2 ’ a ||^ 1/2 + ||(P 3 ’“ + (d“V*77) • u )|| 2 + k 2 ||0“S||§ 


2N 


m) + (p(£°iv}) + tMmu] + £°[ V d + p(£[ti mu] +IF 

3= 0 


4N—2j 


, (6.26) 


where P is a universal positive polynomial with P(0) = 0. Next, we use Cauchy’s inequality, 
Proposition lE.3[ trace theory, and the estimates || J|| L oo + ||A7|| Loo < 1 (which follow from Lemma 
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IE. II and the bounds on £[j?](T)) and ([6. 261) to bound, for any 0 < E\ < 1, 


[ Jd a F 1 ■ d a u — [ 

.hi J e 


• d a u - / d Q F 2 -d a u + I JF l,a ■ d a u — I F 2 ’ a ■ d a u 


1 ,a 


L 


p2,a 


< 


+ [ d a r ] (d a u ■ AT) + [ (nd a E - (F 3 ’“ + • u))0“r/ 

J E is 

Ei ||<9“u|| 2 + 1 (ll^llg + + \\d a F 2 \\l + || F^Wl + \\F 2 ’ a \\ 2 _ 1/2 ) 


+ ||E 3 ’ q + (0°V*77) -u||o + k 2 ||5“H||q < £l ||<9 Q u||? + - H^Hg 

£l 

1 / 27V II 

+ - mi + (p(f 0 M) + e=oM)(£M+«°[ d]) + p(mm + «£ 


3=0 


2 

AN—2j 


(6.27) 


(6.28) 


For the remaining H terms on the right side of (16.181) we estimate 

r ^ r 2N 

/ Kd a ZaA*d a r]<- Ka\A*d a ri \ 2 + — ^ mi 

JYl ^ JYj ^ j_ o 

with the aim of absorbing the A* 3“?7 term onto the left of (16.181) . Next we handle the (E 3 ’" + 
(<9“V* 77 ) • ■u)(ctA*< 9 “? 7 ) terms on the right side of (16.181) . We integrate by parts and then use 
(16.231) of Lemma 16.21 to see that 


2 

4N-2j 


[ -(F 3 ’ a + {d a V*r,)-u)aA*d a ri= [ V*(F 3 >° + (<9“V*r/) • u) ■ aV*d a rj 
is is 

<^||V,(F 3 ’“ + (^V,r7)-«)||o+ [ ^\^r,\ 2 <S°[r,]S[u]+£[ V ]T>[u] + 


/J |V ^ |a - 

(6.29) 


It remains to handle the term 


J (d a V*ri) ■ u(-d a rj + aA*d a rj). (6.30) 

We have that 

— [ (< 9 “V* 7 ?) • ud a r] = - f l-u ■ V* |< 9 q 7/| 2 = j \ div* u\d a r]\ 2 
is is 2 ,7s 2 

< lldiv. u\\ L „ £ 1| 8”>)| 2 < 3M(r) ^ 1| 8”>)| 2 . (6.31) 

Similarly, 



< || V*ii|| Loo | |V*3 a r?| 2 < 3[u\(T) J 1 IV^I 2 . (6.32) 
Combining (16.301) (16.321) . we find that 

[ (3 Q V*r/) • u(-d a r] + aA*d a rj) < 3[it](T) [ \ |<9"r/| 2 + ^ |V*<9“?7| 2 . (6.33) 

7e 7e 2 2 

Now we employ the estimates (16.251) . (16.271) . (16.281) . (16.291) . and (16.331) in (16.181) . By choosing 
£1 small enough (but universal), we may absorb the £1 ||5 Q u || 2 term in (16.271) onto the left. This 
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results in the differential inequality 


d ( f pJ 


dt 


1 


1 2 , O’ 


gt ( / VI d a <+ / ^\d a ^ + -\W^\+C\\d a u\\{ + - / \V*d a rj\ + a \A*d a r]\ 


12 K 


< 

rs_/ 


(i + ||a i (pJ)(pJ)- 1 || i „ + 3H(r» (jf |a“u| 2 + jf i |a“,| 2 +| |v.9”,,| 2 ) 

2N 


+ Mt) + (P(£°[v}) + <£ooM)(£M + £>]) + p(z[t,])t>[u] + «E |F 

l=o 


4N-2j 


(6.34) 


Notice that 


sup || d t (pJ)(pJ) 1 || LOO < 3r[r/](T) + <£\p](T). 

0<t<T 


(6.35) 


We may then apply Gronwall’s inequality to (16.341) and sum over a to deduce the bound 


2 TV 

E 

1=0 


U 


L™H° 


2N 


+ 


a)v“- 2 A 


L^H 1 


+ 


dj/n 


L oo H 4N - 2 3 


+ k E n 77 


3=0 


L^H 4N - 2 j 


+ cr 


<9^ v* 7 ? 


L^H 4N ~ 2 i 


+ cr <9^V* 7 7 

2 TV 


L™H 4N ~ 2 i 


+ e 


7 t 


f Mt) + (P(£°M) + «coW)(f[»] + £>]) + P(£M)P[u] + K £ ||a| 
■'» V j-o 


2 

4IV-21 


dt, 

(6.36) 


where 7 = (7 (1 + 3[r/](T) + (£[p](T) +J[u](T))). Then (16.241) follows from (16.361) and (15.50 . □ 

Next we employ various elliptic estimates to gain control of all derivatives of ( u,r 7). 
Proposition 6.4. Assume (16.50 - (16.120 . Then for 0 < T < T* we have the estimate 


2 N 


2 TV 


£ K 

3=0 


u 


+ \\u r 2 ttAN + 




l 2 h 2 

T 3=1 


L 2 t H 4N ~ 2 1 +1 


+ 


pJc) t 2Ar+1 « 


l 2 t { 0 h 4 y 


2N 

E||^ 

1=0 


L9pH 4N ~ 2 i 


+EH 

2N 

E^ 1 5 i v *7 


1=0 


L o°h 4N - 2 1 


2N+1 


+ °’ 2 IMIl 2 /74JV+3/2 + 11 <9*711 L2 J/4AT^l/2 + £ ^7? 


1=2 


L 2 t H 4N ~ 2 3+ 2 


2 TV 

+k E 

1=0 


L 2 T H 4N ~ 2 i 


&j-A *T) 


L 2 T H 4N ~ 2 i 


+ T 


<e lT 4"h,')0] + sup T'oo(i) 

0<t<T 

sup £°[77(f)] + e 7T f 1 (j- 2 (t) + (P(£: 0 [7;]) + P(C:oo[p]))(^M+£: 0 [7]) + P(T[r;])^M)^ 
<t<T Jo 

+ [ £[u\— s a [r]\dt + K 2 3, (6.37) 
lo o' 

where P is a polynomial with positive universal coefficients such that P( 0) = 0 and 7 = 

c(i + a[i7](r) + (SW(T)+a[u](r))). 

Proof. The argument is very similar to one used in ns, so we will only provide a sketch. Let 
us write 2 to denote a term of the same form as the right hand side of (16.240 . From line to line 
we will let the polynomials vary as well as the universal constant C > 0 appearing in 7. 
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First we use trace estimates to bound 
2 N 

£ IN“ 


j=0 


2 A 


L |^4iV-2 3 '+l/2( E ) 


< e 


3=0 


L 2 T m 


< z. 


(6.38) 


Then we use elliptic regularity for the Lame system with Dirichlet boundary conditions, Propo¬ 
sition IE.41 to deduce that 


2 N 


\ u \\ 2 L 2 T H iN + E In' 

3 = 1 


It 


L 2 H 4JV-2j + 1 


< z. 


(6.39) 


Notice here that we have only used the estimate (IE.61) in order to avoid the introduction of the 
term |N|l 4 iv+i/ 2 - This causes no difficulty in getting 41V — 2j + 1 estimates when j = 1,..., 2 N 
but prevents us from obtaining an estimate of | N11 4 jv +1 at this point. Instead we sum the elliptic 
estimate for d a u obtained from (16.131) with a € N 2 and |a| < 41V — 1 to obtain the estimate 

2 




< Z. 


(6.40) 


Now we derive the L°° in time estimate for u and its time derivatives. Lemma ID. II provides 
us with L°° in time estimates at lower regularity: 


2IV—1 

£ 

3 =1 


din 


L <x, H iN-2j 


2N 

S«oN + £||a?' 

3 = 1 


LlH iN ~ 2J+1 ’ 


(6.41) 


and hence (16.391) and the bound of j|<9f A N||Ef/o provided by (I6.24P imply that 


2 N 


E || d t i 

3 =1 

Similarly, Lemma ID. II and (16.381) imply that 


Ltj?H iN - 2 i 


< Z. 


(6.42) 


i 

lL 5 ?H 4 Ar - 1 / 2 (S) ~ 4 No] + ^2 || 

3=0 


u 


L 2^4iV-2,+l/2( S ) 


< z. 


We may then use (16.431) and the elliptic estimate of Proposition IE. 41 to bound 


(6.43) 


\ u \\l™H an ~ \\pdtu\\ L cc,H 4N 


-2 + \\F 


\l?^h^ n ~ 2 T INI 


L oo H 4JV- 1 /2( E ) 


&\p\(T) \\ d M\ L ^H ^-2 + sup Fooft) + |NIIl 2 ?H 4JV - 1 / 2 (S) 

T 0<t<T T v ; 


(l + £[p](T))Z + sup .Foo(t)<2 + sup ^(f). (6.44) 

0 <t<T 0 <t<T 


Note that in the last inequality we have used the bound <£\p](T)Z < Z, which holds since we 
may increase the constant C > 0 appearing in 7 . 

Next we derive the L 2 in time estimates for 77 and its derivatives. We first use the dynamic 
boundary condition and (I6.40p to estimate 


a 


lL2^4JV+3/2 


< —2 II — „ 2 


a 


l 2 t ho + II ^* 1 i\\l 2 t lrw-1/2 


< 


T sup f°[r i(t)\ + Z+ [ ((l + P(£°[v])\\V*v\\lN-i/ 2 £ [ u ]) dt 

0 <t<T .In \ * 


< 


0 <t<T 


T sup £°[r)(t)\ + Z + / ( (1 + P(£°[rj\)— £ a [rj\£[u] ] dt. (6.45) 


a 
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Then we use the parabolic modification of the kinematic boundary condition to get improved 
L 2 in time estimates for dfr] for j > 1. Standard heat equation estimates yield the bounds 




2N+1 

+ | 5 * ?? 

3 =2 


L 2 t H 4JV-2J+2 


<Z + [ £[u\-£ a [rj\dt + k 2 3- (6.46) 

Jo a 


We then combine the estimates (16.241) of Proposition 16.31 with (16.391) . (16.401) . (16.421) . (|6.44l) . 


(I6.45p . and (16.461) to deduce all estimates in (16.371) except that of 

recover this estimate we simply appeal to (16.131) with d = d\ 

2 


then allows us to estimate 


pJd 2N+1 u 


this yields (16.371) in the form written. 


l 2 t ( 0 h^Y 


Jd 2N+1 u . To 

L 2 T ( 0 my 

,2N A standard duality argument 
in terms of all of the existing terms in (16.371) : 

□ 


Finally, we combine Propositions 16.31 and 16.41 in order to obtain n —independent estimates. 


Theorem 6.5. Assume (16.51) (16.121) . Assume also that 0 < k < min{l, cr+, <r_}. There exists a 
universal constant 0 < 5\ < 5 q/ 2 (where <5o is from Proposition \E.4\ > and a 0 < T\ = T\(£q , a) < 
1 such that if 0 < S < 5i in (J6.5D (16.121) and 0 < T < min{Ti,T*}, then 


sup {£[u{t)\ +£ a [r]{t)}) + [ (V[u{t)\+ pJd 2N+l u(t) +V a [rj(t)])dt 

<t<T Jo (o H 1 )* 


0 <t<T 
2N rT 

+*E/ 


^fV*7 l(t) 


4N—2j 


+CF 


dt 


dt < P(£q [n 0 ,??o])+ sup JFoo{t)+ JF 2 {t)dt 
4N—2j 0 <t<T Jo 

(6.47) 


for a positive universal polynomial P such that P( 0) = 0. 

Proof. We first notice that 

7 = C (1 + J[ V ](T) + £[p](T) +Z[u](T))) <C( 1 + P(4")), (6-48) 

and so we can make T 2 small in terms of £q , cr, and a universal constant so that e^ T < 2, 
VTP((B[p\(T)) < 1, and \[T < a whenever T < min{T 2 ,T*}. Then Proposition 16.41 yields the 
estimate 


sup {£[u(t)} +S a [p{t)}) + 
0 <t<T 


{V[u(t)} + pJdj N+1 u{t) 


+ V a [r](t)})dt 

(oH 1 )* 



2 

+ <7 

4N-2j 


di A * v(t) 


2 

4N—2j 


dt < C£(f [u 0 ,r]o] + C sup Poo(t) 

o <t<T 


+ C [ P 2 {t)dt + Cn 2 f, + VTP 

Jo 


sup (£[u(t)\ +£ a [p(t )] 
0 <t<T 


+ C sup £[??(£)] 

0 <t<T 


V[u(t)\dt 

(6.49) 


for every T < min{T 2 , T*}, where P is a universal positive polynomial such that -P(O) = 0, and 
C > 1 is a universal constant. 

Since 


sup 2(q(t)] < 5 < Si, 
o <t<T 


(6.50) 
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we may choose a universal <5i > 0 such that C8\ = 1/2. We may then absorb the last term on 
the right side of (|6.49l) onto the left. This yields the estimate 


sup (£[u(t)] +£ a [v(t)}) + \ [ (V[u(t)]+ pJd 2N+1 u{t) 

0<t<T 2 Jo 


(oH 1 )■ 


+ V a [rj(t)])dt 


2 TV ..t 


+ k Y 2 f ,9 t v * r ?( t ) 

j=o Jo 


AN—2j 


+ C 7 


di A *v(t) 


AN—2j 


dt 


< C8o[u 0 ,rio] + C sup .Foo (t) + C F 2 (t)dt + Ck 2 3 
0 <t<T Jo 

+ VfP I sup (£[u(t)\ + £ a [p(t)}) ) (6.51) 

\0<t<T J 

for every 0 < T < min-fT^T*}. 

We may view (16.5111 abstractly as an inequality of the form 

X(T) < CZ(T) + VTp(X(T)) (6.52) 

for X, Z : [0, min{T 2 , T*}] —>• [0, oo) continuous functions such that Z is non-decreasing and 
X(0) < CZ{ 0). By continuity we know that either X(T) < 2 CZ{T) for all 0 < T < min{T 2 , T*}, 
or else there exists a first time 0 < T 3 < min{T 2 ,T*} such that X(T 3 ) = 2CZ(T 3 ). Plugging 
this equality into (16.521) implies that 

2 CZ(T 3 ) = X(T 3 ) < CZ(T 3 ) + y/T 3 P(X(T 3 )) = CZ(T 3 ) + v^P(2CZ(T 3 )), (6.53) 

and hence 

CZ(T 3 ) < y/ ¥ 3 P(2CZ(T 3 )). (6.54) 

From this we deduce that 




nun 


where 


Zmax — 


2 P{2CZ(T 3 )) 2 ze[o,2Cz max ] P(z ) ’ 


max Z(T) = Z(mm{T 2 , T*}) < P{£%) + C. 

0<T<min{T 2 ,T,} 


(6.55) 


(6.56) 


The last estimate follows because of assumption (16.111) . from which the universal constant C > 0 
comes, and (16.4p combined with the bound k < min{ 1 , <t_|_ , <r_}. Since P( 0) = 0, we then hnd 
that 

\fp 3 > ^ min _ -X—- =: > 0. 


min —— 

2 ze[0,2C(P(£g)+C)} 

This leads us to define 

T 1 = T 1 (£Z,a) = mm{T 2 ,T 4 } 

so that if 0 < T < min{Ti, T*} we have the estimate 

X(T) < 2CZ(T). 

Removing the abstraction, this implies that 

sup {£[u{t)\+£ a [rj{t)}) + ]- I (V[u{t)} + pJd? N+1 u(t) + V a [r](t)})dt 

0 <t<T 2 Jo ( 0 H 1 )* 


(6.57) 

(6.58) 

(6.59) 


2 N rT 2 
3t’v* ri(t) 


+ «E/ 

i=o Jo 


AN—2j 


+ a 




AN—2j 


dt 


< 2C£q [uo,rjo] + 2C sup F 00 (t) + 2C / P 2 (t)dt + 2C'k 2 3 (6.60) 
0 <t<T Jo 


for 0 < T < min{Ti,T*}. This and (16.41) yield (|6.47[) . 


□ 
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7. The two-phase free boundary Lame problem 


Our goal in this section is to produce a solution to 


" pdtu — div _4 §,477 = F 1 

in Q 

dtp = u ■ M 

on £ 

—S_AuJ\f = —a + A*r]J\f + F+ 

on £_|_ 

— A f = cr_ A*r/Af — Fj 

on £_ 

N =0 

on £_ 

u- = 0 

on 

M',0) = u 0 ,p(-,0) = 770 



by passing to the limit n —> 0 in the k —approximation problem (15.11) . 

Our strategy is as follows. First, we will combine the local existence result of Theorem 
15.11 with the a priori estimates of Theorem 16.51 to produce solutions on a time interval that is 
independent of k and that satisfy k —independent estimates of the form (16.471) . Second, we will 
pass to the limit k —>• 0 to recover a solution to (ED- 


7.1. Local existence on At—independent intervals. Our goal now is to combine Theorems 
15.11 and 16.51 We first describe some assumptions on the data and forcing terms that will be 
needed. 

We say the data and forcing satisfy on the time interval [0, T*] if the following hold. 


The initial data dju(-, 0) and 0) for j = 0,..., 2N are given as before and (7-2) 
satisfy the compatibility conditions (15.61) . 

The data satisfy £[ 770 ] < 5/2, where £[ 770 ] is given by (1A.10D . (7.3) 

f T * 

The forcing terms satisfy sup Foo{t) + / T 2 (t)dt < P(£q ) 

o<t<T* Jo 

for a universal positive polynomial P such that P( 0) = 0. (7-4) 

The function p is defined on [0, T*] and satisfies (fT 2 l) and flO}. (7.5) 


Theorem 7.1. Assume that 0 < n < min{ 1, cr+ , cr_} . Let 0 < 5i he the universal constant 
and 0 < T\ = Ti(£q ,< 7 ) be from Theorem 1 6. 5\ Assume that the data and forcing satisfy ip(5i) 
on the time interval [0, T*]. Then there exists 0 < T 2 = T 2 (£q,(t) < T\{£q,ct) such that if 
0 <T < min{T*,T 2 }, then the following hold. 

( 1 ) A unique solution (u, 77 ) to (|5.1I) exists on the interval [0, T] and achieves the initial 
data. 

(2) The solution satisfies the estimate 

21[r](T) +©[t/](T) < 00 , (7.6) 

where 21 [r] and 23[ 77 ] are as defined in (15.91) . 

(3) We have the estimate 


sup (£[u{t)\+£ a 
0 <t<T 


fo(*)])+ [ (£[«(*)]+ pJd? N+ 1 u{t) 

Jo 


+ V a [ 77 (f)] )dt 
(0 hw L,WJ; 



2 

+ C7 

AN—2j 


djA* 7?W 


4N—2j 


dt < P(£q [uo,7?o])+ sup Too(t)+ / P 2 (t)dt 


0 <t<T 


(7.7) 


for a universal positive polynomial P such that P( 0) = 0. Here we recall the notation 
V defined in (16.11) . 

(4) We also have the estimates 


£[ 77 ] (T) < St 


(7.8) 
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and 

3[r,](T)+3[u}(T)<P(£Z) + 1 (7.9) 

for a universal positive polynomial P such that P( 0) = 0, where £ is defined by (IA.9D 
and 3 and J are defined by (15.101) . 

(5) One of the following is true. Either T 2 = T\, or else there exists a universal positive 
polynomial P satisfying P( 0) = 0 and a universal constant C > 0 such that 


o5i 

C(1 + P(£S )) 


< T 2 . 


(7.10) 


Proof. For r > 0 let (5(r) denote the proposition that the following three conditions hold. First, 
a unique solution to (15.11) exists on [0, r] and achieves the initial data. Second, the solution 
satisfies (17.6p with T replaced by r. Third, the solution satisfies (17.81) and (17.91) with T replaced 
by r. Define the set 

91 = {r € [0, rnin{T*. T\ }] | 6(r) holds}. (7.11) 

Theorem 16.51 guarantees that di < do/2, where <5o > 0 is defined by Proposition IE. 41 We may 
then apply Theorem 15.11 to see that 91/0. Let 7/ = sup91 € (0, T*]. 

If 7/ = min{T*,Ti} then we set T 2 = Tf, and we are done. Indeed, the hypotheses of 
Theorem 16.51 are satisfied, and so (17.71) follows. We may assume then that T/ < min{T* ,T\} 
throughout the rest of the proof. 

If 7/ = max 91, then a standard continuation argument, employing Theorem 15.11 to extend 
the solutions, yields a contradiction, and so we deduce that 7/ ^ 91. This means that ©(71 yf) 
fails. We claim that in fact it is only the third condition in £5(7/) that can fail, i.e. the first 
two conditions remain true at 7/. 

We know that £5(7/ — e) is true for e sufficiently small. The a priori estimates of Theorem 
16.51 are then valid and provide e— independent control of £(7/ — e) in terms of £ff. We may then 
use Theorem 15.11 to extend the solutions to T 3 = 7/ — e + To( k, a, 5, £/f). When e is sufficiently 
small we have that T 3 > 7/, and so the first condition of ©(7/) must hold. Additionally, the 
functional setting of Theorem 15.11 guarantees that the second condition of ©(7/) must also hold. 
We deduce then that it is the third condition of ©(7/) that fails at 7/, proving the claim. 

From Lemma lD.il Theorem 16.51 and (17.41) we may estimate 

3fr](?/)+3[u](7Z) < £0 [u 0 ,Vo]+ [ (£[u(t)]+£°[r)(t)])dt<£$[u 0 ,r) o]+/ P(£o)dt, (7.12) 

Jo Jo 

and since 7/ < 7\ < 1, we then have that 3[rj\(T 2 ) + 3[u](T 2 ) < P(£ q) + 1. This means that 
estimate (17.91) remains true at 7/, so it is actually estimate (17.81) that fails at time 7/. Arguing 
similarly, we deduce that 


Si = -CM?/)] 


MTn)\\l N _ 1/2 <2[r,o\ + C 


AN + \\dtri{t)\\l N -i )dt 


<h + 1 

~ 2 a 


[ T ” £ a lv(t)]dt < | + ^(1 + P(£%)), (7.13) 
Jo la 


and hence that 

cr5\ 

--- < Tat 

2C(1 + P(£ 0 ff )) - 

To conclude the proof we then set T 2 = Tr^/2, apply Theorem 16.51 and use (I7.14|) 

(fTroD . 

7.2. Sending k —>■ 0. Our aim now is to use Theorem 17.11 to send k —>■ 0 in 
produce a solution to dm 


(7.14) 

to produce 
□ 

in order to 


Theorem 7.2. Let 5\ > 0 be the universal constant from Theorem 1 6. 71 and assume that the 
data and forcing satisfy fp(5i) on the time interval [0, 7 1 *]. There exists a 0 < T 3 = Ts(£q) such 
that if 0 < T < min{7*, Ti}, then the following hold. 
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(1) A unique solution ( u,r 7 ) to (17.11) exists on the interval [0, T] and achieves the initial 
data. 

( 2 ) We have the estimate 


sup (£[u(t)}+£ a [rj(t)])+ 

0 <t<T 


(V[u(t)}+ pJd 2 t N+1 u(t) 


(off 1 )’ 


+ V a [ri(t)})dt 


< 


P(£q [t/ 0 ,7?o]) + sup -7 r oo (0+ / F 2 (t)dt 
0 <t<T Jo 


(7.15) 


for a universal positive polynomial P such that P( 0) = 0. 

(3) We have the estimates 


mi?) < Si and 3[rj\(T) + 3[u\(T) < P(£%) + 1 (7.16) 

for a universal positive polynomial P such that P( 0) = 0, where £[ 77 ] is defined by (IA.9D 
and 3[rj\ and 3[it] are defined by (15.91) . 

(4) We have the estimate 




2N+1 

\In-1/2+ Y I ^(t) 

3=2 


4N—2j+5/2 


dt 


< 


P{£{ 


r ) + sup JFooft) + f T 2 (t)dt 
0 <t<T Jo 


(7.17) 


for a universal positive polynomial P such that P( 0) = 0. 

(5) We also have the estimate 


o<t<T ^^^iv+ 1/2 < ex P (cT f lkWII^4iv+i/2 (s) dt 


x ( ll^oII 4 JV+ 1/2 + T J o \\ u ( t )\\H iN + 1 / 2 CE) dt) ■ (7.18) 


Proof. We divide the proof into two steps. In the first we will initially prove the theorem with 
the weaker condition that T 3 = Tfi £^, er), i.e. that T 3 is allowed to depend on a as well as the 
data, but not on k. In the second we will use the first step and some auxiliary arguments to 
remove the a dependence. 

Step 1 - T 3 = Tfifif.o). 

For each 0 < k < rninjl, <r + , <r__} Theorem [T7T] provides us with a pair (u K ,r) K ) solving (15.11) on 
(0, T 2 ) and satisfying the conclusions of the theorem. We shall consider k to index a sequence of 
values chosen in (0,min{l, cr_|_, <r_}) that decrease to 0. The k —independent estimates of (17.71) 
show that the sequence {(u K , r] K )} K is bounded uniformly in the function space determined by the 
first line of (Ell) (i.e. the left hand side of the inequality with k = 0). These uniform bounds 
allow us to argue as in Theorem 6.3 of m, using weak and weak-* compactness arguments 
together with interpolation and lower semi-continuity arguments, to extract a subsequence (still 
denoted by k) such that 

diu K -A d}u in C°([0, T\:H 4N -^-' 2 (n)) for j = 0,1, 2 

tfrh -> ^ in C°([0, T\. (£)) for j = 0,1, 2, 

where (u, rj) satisfies (17.71) and achieves the initial data. Note here that the convergence (17.191) 
can be improved to a larger range of j and to higher regularity for the various 77 terms. We state 
(17.191) as is because it is more than sufficient to pass to the limit in (15.11) and deduce that (u, 77 ) 
are a strong solution to (17.11) . The estimates (|7.16l) follow from (|7.8D and (17.91) using similar 
lower semi-continuity arguments. 

It remains to derive the improved estimates. This entails improving (17.71) to (17.151) by obtain¬ 
ing an estimate for Hulljy ^ 4 jv+i and also proving (17.171) . To accomplish the first we will need 
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(17.1811 . This estimate is proved in Step 1 of Theorem 5.4 of [TTJ by using estimates developed 
in [7] for solutions to the kinematic transport equation (I7.25H . 

With (17.1811 in hand, we use (17.711 and (17.411 in conjunction with the simple estimate 


I u IIh 4jv + i /2(s) 




< 

H 3/2(S) ~ 




in order to bound 


SUP ||»7(*)|l4JV+l/2 ~ eX P 
0 <t<T 


<CTP (£° 0 


£S + TP(££ 


Then if T 3 < T 2 is taken sufficiently small, we may bound 

sup h(t)\\l N+1/2 < P{£o). 

0 <t<T 

We then appeal to the elliptic estimate (IE.71) of Proposition IE. 41 to bound 


(7.20) 

(7.21) 

(7.22) 


j-T rT 

J \Ht)\\l N+1 dt < J (||p9tu(i)llLv-i + ||E 1 (t)||4 JV _ 1 + ||u(t)||^4iv+i/2 {E) ) dt 

pT 

+ 0 ^ T hml N+ y, l (\\P d Mt)\\l + ||'^ ll (^)||2 + IN*)lltf7/2( E) ) dt. (7.23) 
Hence Q>, (ED, and (HZD allow us to bound 

[ \\ u {t)\\lN+i dt ~ [ {D[u(t)]+P 2 {t)) dt + T(l + P(£g)) sup (£[u(t)\ + Pooit)) 

Jo Jo 0 <t<T 

<P(Sq) + sup T-oo(t)+ [ T P 2 (t)dt (7.24) 
0 <t<T Jo 

once we further restrict T 3 so that T^l + P(£^)) < 1. Summing (17.2411 and (17.711 then yields 
(I7T5D . 

Finally, we prove the improved estimates (I7.17H . For this we will use the fact that (u,p) 
satisfy the kinematic equation 

dtp = U 3 + V*?] • u on E. (7.25) 

We may use this equality as in Theorem 5.4 of EH, using the usual estimates of products in 
Sobolev spaces, to deduce that 


2 N 

sup y Mvit) 

0<t<T “i 11 


4Af-2j+3/2 


+ 


27V+1 

't''/ll4V-l/2 + X] || 
3 =2 


4Af-2j+5/2 


dt 


<P sup (£[u(t)] +£ a [p{t)]) + / (V[u{t)\+V a [p(t)])dt\ (7.26) 

\ 0 <t<T Jo J 


for some universal positive polynomial with P(0) = 0. Then from (17.2611 . (I7.15H . and (17.411 we 
deduce that (17.171) holds. 

Step 2 - Improvement to T 3 = T^(£q) 

With the Theorem in hand for a T 3 = T^{£q,g), we can employ a continuation argument in 
order to remove the dependence on a. The argument is similar to the one used in Theorem 17. II 
so in the interest of brevity we will only point out the key point. This lies in the fact that our 
result from Step 1 requires fp(di) to hold on [0,T*], which in turn demands that ||^ 0 II 47 V— 1/2 < 
5i/2. The estimate (17.1711 provides us with an estimate of sup 0<t<T \\dtp(t)\\‘l N _ 1 ^ 2 , which when 
coupled to the fundamental theorem of calculus, allows us to estimate £[ 77 ](T). Using this, we 
can prove that the theorem remains true on a time interval [0, min{T*, T 3 }], where either T 3 = T* 
or else T 3 is when \\p(T 3)\\‘^ N _ 1 ^ 2 > We may then argue as in ( 1 7. 1 4 1 ) . using the estimate for 
\\dtp\\1 N _i / 2 in place of Lemma fD.ll to show that T 3 is bounded below by a positive function 
of £q that is independent of a. Hence T 3 = T^(£q). □ 
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Remark 7.3. The improved g estimates of (17.1711 are the key to eliminating the dependence of 
the existence interval on a. These in turn depend on the structure of the equation dtq = u-J\f in 
place of the n approximation. We thus see the importance of passing to the limit k —>• 0 before 
attempting to remove the dependence on a. 


8. The transport problem 
Our goal in this section is to study the transport problem 

f d t q + ( A T u - Kd t 9e 3 ) • Vg + div A (u)q = f 
= g 0 

where u, r) (and hence A, etc), and / are given. To simplify 
initially study the more general problem 

d t q + v - Vq + cq = f, 


inflx(0 ,T) 
in 12, 

the structure of m we will 


( 8 . 2 ) 


where v, c, and / are given. 

The key to the analysis of a transport problem (18.21) on a domain with boundaries is the 
behavior of v at the boundary. For the v arising in ([8.111 with (u, rj) satisfying (17.111 we have the 
following crucial identities: on we have 

v ■ e 3 = (A t u — KdtOe 3 ) ■ e 3 = K(u ■ JAe 3 — dtq) = K[u± ■ J\f± — dtr]± ) = 0, (8-3) 

while on £& we have 

v ■ 62, = ( A t u - Kd t 9e 3 ) ■ e 3 = K(u- ■ JAe 3 - bid t fj + - b 2 dtfji) = 0 (8.4) 

since b\ = b 2 = 0 and u- = 0 on £&. Note that the condition u- = 0 on E& could be weakened 
to u 3 - = 0 on S;, since JAe 3 = e 3 on £&. The upshot of these identities is that we may assume 
that the vector field v satisfies 

v ■ e 3 = 0 on <912. (8-5) 

Notice that we do not couple the transport equation to boundary conditions, and in fact the 
equations in 12+ and 12_ are decoupled from one another. This allows us to solve the equations 
in each domain separately. To this end we will let T = 12 + or T = 12_ and discuss the transport 
equation in T. Let us now define various functionals that will appear in our analysis of the 
transport problem. 

For a function g we define 


and 


He [g] 



2 

4N-2j+l 


and %K e \g] 


2N—1 



2 

AN-2j-\ ' 


( 8 . 6 ) 


o. d [g\ 


2 

AN—l 


2N+1 

+ £ 

3= 2 


dig 


2 

4N-2J+2 


and <R d [g\ 


LiV-i 


+ J 2 In 9 


3= 1 


2 

AN—2j ' 


(8.7) 


We will assume that v, c, and / satisfy 

sup (JR e [v(t)] +9\ e [c(t)] + 9t e [/(£)]) + f (Rd[v(t)]+Rd[c(t)]+ftd[f(t)])dt < 00 . (8.8) 

0 <t<T Jo 


8.1. Solution by characteristics. Consider the transport problem 

I d t q + v ■ Vg + eg = / in T x (0, T) 

\<?(-,0) = g 0 in T. 

Here we assume that v, c, and / satisfy (18.811 . In particular the usual Sobolev embeddings 
require that v, c , and / are all C 1 (T x [0, T]). 
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To produce a solution to (18.91) we use the method of characteristics. We let Ct(x) denote the 
solution to the ODE 


d t (t(x) = v(Ct(x),t) 
Co(x) = x 


( 8 . 10 ) 


for x £ T and t € [0, T\. The identity (18.51) is essential here, since it guarantees that for each 
t £ [0, T], ■ T —» T is a diffeomorphism. Let us define the map uj : [0, T] 2 xTaT via 


Then the solution to (18.91) is 


u(s,t,x) = Cs(C t V))- 


( 8 . 11 ) 


q(x,t) = q 0 (cu(0,t,x))exp 


/ c(cu(s, t, x), s)d. 

Jo 

+ / f(u(s,t,x),s) 
■Jo 


exp 


/ c(u!(r,t,x),r)dr 
J S 


ds. (8.12) 


While this explicit form of the solution is nice, it is not convenient for making higher regularity 
estimates. In particular it is not immediately obvious from (18.121) that q belongs to the space 
defined by (1A.2D and (|A.5[) . and it is not clear that we can justify applying d a to (j8.9D and 
performing a priori estimates. The usual solution to this difficulty is the Friedrichs mollification 
method: we first solve a mollified version of (18.91) so that the solution is smooth enough to 
justify the a priori estimates, then we derive the a priori estimates in a manner independent 
of the mollification parameter, and then finally we pass to the limit. This works well when 
ST = 0, but the mollification procedure runs into technical obstructions when ST ^ 0, which 
is the case here. 

Two options then present themselves. The first is to modify the mollification procedure in a 
manner that makes sense in our T but does not destroy the structure of the a priori estimates. 
The second is to transfer the problem (18.91) to a new problem on a set without boundary in 
such a way that the estimates from Friedrichs’ method carry over to (|8.9D . We have chosen to 
go for the second option since Friedrichs’ method is so well-known. Our goal then is to justify 
the transfer of the problem and the estimates. 


8.2. Transfer. Consider the transport problem (18.91) where we only assume for now that q$ £ 
L 2 (T) and / £ L 2 ([0, T]\ L 2 (T)). We say that q £ L 2 ([0,T];L 2 (T)) is a weak solution to (18.91) 
if 



-q (d t p + div(n^) 


cp) 




Qop('iO) 


(8.13) 


for all p £ C^(T x [0, T)). The identity (|8.13l) is clearly satisfied by any regular solution to 
(18.91) . and in particular by the solution given by (|8.12l) . 

We have the following simple lemma on the uniqueness of weak solutions, which we state 
without proof. 


Lemma 8.1. The following are equivalent. 

(1) For every f € L 2 ([0, T]\ L 2 (T)) and qo £ L 2 (T), there exists at most one function 
q€ L 2 ([0,T];L 2 (T)) that is a weak solution to (18.91) . 

(2) If q £ L 2 ([0, T]; L 2 (T)) satisfies 



-q ( d t p + div(v^) 


cp ) = 0 


(8.14) 


for every p £ C x c (T x [0,T)), then q = 0. 


Our next goal is to verify that the second item of Lemma 18.11 holds, which means that weak 
solutions to (18.91) are unique. To this end we first study the adjoint problem determined by 
what appears in parentheses in ([8.141) . 
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Let ip G C %°(T x (0,T)). We want to find p G Cp (T x [0, T)) satisfying the adjoint problem 

I d t ip + div^) — c<p = ip in T x (0, T) . 

U(-,T)=0 inT. 


Note that this is a terminal value problem; we seek to solve this so that <p can be taken to be 
compactly supported in T x [0,T). 

To solve the adjoint problem (18.1511 we use Q given by (18.1011 to reduce to an ODE along the 
characteristics. Recall the function u : [0,T] 2 xT->T given by (18.1111 : we may use it to derive 
the solution 


<p(x,t) = -J ip(u(s,t,x),s) 


exp 


(divu(u i(r,t,x),r) — c{ui(r,t,x),r )) dr 


ds. (8.16) 


From this formula and the inclusion ip G Cp°(T x (0, T)) it is easy to see that <p G Cp (T x [0, T)), 
as desired. We have thus proved the following lemma. 


Lemma 8.2. Let ip G C £°(T x (0, T)). Then there exists p G C] (T x [0, T)) solving (18.15H . 


With Lemma 18.21 in hand we can prove that weak solutions to (18.91) are unique. 

Proposition 8.3. Let f G L 2 ([0, T]; L 2 (Y)) and qo € L 2 (Y). Then there exists at most one 
function q G L 2 ([0, T]\ L 2 (T)) that is a weak solution to (18.911 in the sense of (18.131) . 


Proof. To prove uniqueness we will show that the second item of Lemma fe. II holds. Suppose that 
q G L 2 ([0,T];L 2 (T)) satisfies (|8.14l) for every p G C \{T x [0,T)). For any ip G C %°(T x (0, T)) 
we may use Lemma [8721 to find <p G Cp(T x [0,T)) solving (18.151) . Using this p in (18.1411 yields 
the equality 

/ /# = 0. (8.17) 

Jo Jr 

Since ip was arbitrary, we deduce that <7 = 0. Hence the second item of Lemma 18. II holds. □ 


Next we define the extended problem, which is easier to handle with Friedrichs’ method. Let 
T = T 2 x I denote the extended domain in which we will pose the extended problem. Let E 
denote a Sobolev extension operator such that E : H m ( T) —>• H m (T) for every m = 0,... ,M, 
where M > 0 is large enough to surpass every regularity index in the energy and dissipation 
defined in (1A.2H and (1A.5I) . We then define 

v = Ev, c = Ec , go = Eqo, and / = Ef. 

We say that q G L 2 ([0, T]; L 2 (T)) is a weak solution to 

id t q + v ■ S7q + cq = f in T x (0,T) 

\?(-,0) = g 0 inT 

if 

[ -q (d t p + div(vp) - cp) = [ [ pf + [ q o p(-,0 ) (8.20) 

r Jo Jr Jr 

for all p G Cl (T x [0,T)). 

Lemma 8.4. Suppose that q G L 2 ([0, T]; L 2 (T)) is a weak solution to (18.191) in the sense of 
(18.2011 . Let q denote the restriction of q to T x (0, T). Then q is the unique weak solution to 
(18.91) in the sense of (18.131) . 

Proof. Since (18.2011 holds for all p G C 2 (r x [0,T)) it must also hold for all p G Cp (T x [0,T)). 
For such p the equality (18.2011 is identical to (18.1311 because v, c,q o, and / are extensions of v, c, 
qo, and / to T. Hence q is a weak solution to to (18.911 in the sense of (18.13(1 . Uniqueness follows 
from Proposition 18.31 □ 

The upshot of Lemma 18.41 is that if we can produce a solution to the extended problem 
(18.1911 that obeys the estimates we seek for the original problem (18.91) . then we know that those 
estimates are also valid for the solution (|8.12l) given by characteristics on T. This leads us to 
study the extended problem (18.1911 . 



(8.18) 

(8.19) 
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8.3. The extended problem (18.191) . 

8.3.1. The mollified problem. We begin by defining space and time mollification operators. Since 
the horizontal directions in T are periodic, it is convenient to decompose our spatial modifiers 
into horizontal and vertical parts. Let d € C'“(M) be a standard modifier. For n £ N and 
i = 1, 2, let € C'°°(27rLjT) denote the Fejer kernel. Then we define the spatial modifier via 

Ke p g(x) = j^ g{x - y)F{ 1/E \ (yi)F[i/e\ (l/2)^0 (y) dy, (8.21) 

where [z\ denotes the integer part of z. To define the temporal modification operator we must 
first define a temporal extension. For a function g defined on T x (0, T) we first extend to g 
defined on T x R via 


0 


if t <E [0, T] 
if t ^ [0, T\. 

(8.22) 


(8.23) 


g(x,t) = 

Then the temporal modification is 

Kl e g{x,t) = 

The operators K £ p and K £ e satisfy ad the usual properties of a modification operators. 

The mollified problem studied in Friedrichs’ method is 

f d t q E + K! P [(Kfiv) • V(K s £ p q £ )\ + (. Kfc)q £ = Kfif in F x (0, oo) 

\^(-, 0 )=g 0 in T. 

Notice that Kfiv, Kfic and Kfif belong to C°°(M; H iN (T)). Then the theory of linear ODEs 
in Banach spaces provides us with a unique solution q £ € C°°([0, oo); H' iN (T)) to (18.241) . 

Our next goal is to produce e—independent estimates to the solutions to (18.241) . To this end 
we note that for a € N 3 we may apply d a to (18.241) to find that d a q £ solves 


(8.24) 


d t d a q £ + K £ p [(Kfiv) • \/(K sp d a q e )] + (. K^c)d a q e = ff in T x (0, oo) 
d a q £ (-,0) = d a q 0 in T, 


(8.25) 


where 


fa = d a K te-f _ £ ^ ^[(K^v) • V(iFf d a ~^q £ )} + (K^ c)d a ^ q £ ) . (8.26) 

by parts 
(8.27) 


0<f3<a 

8.3.2. Estimates. By multiplying the first equation in (18.251) by d a q £ and integrating 
over r we may derive the basic energy identity 


\d a q e \ 2 + J r ( K l £ d) \d Q Qe\ 2 ~ {K^dwv) \K sp d a q e \ 


= / f £ d a q £ . 


We may use standard Sobolev embeddings and properties of modifiers to estimate 

f £ d a q e < ||/|| 4iV MeWm + 

Then 


Liv + ll c ll 4 iv) II^ILiv • 


d 


1 - 112 


dt imUN 


< 


4 N 


147V + ( 

which leads us to the fundamental estimate 

sup ||9 e (*)|| 4 JV ^ ex P ( C [ (II^WIUlV + MVLn) dt 


+ 11 ^ 114 JV ) lkell 4 V ) 


14 TV 


dt 


90IU1V+ [ ||/( 0 | 

0 <t<T \ Jo J \ JO 

Here the constants on the right-hand side do not depend on e. 

In order to record other estimates we recall the functionals £2 e , given by (18.61) 
9^ given by (18.71) . Here the norms are understood to be computed over T. With the 
(18.301) in hand we may use the equation (18.241) to directly estimate ||II 4 A 7—1 ■ 




i< 


I/II41V-1 + 


-112 


I4AT-1 


+ C 


-H 2 


I4IV-1 


1 - 11 2 

I II4A7 ■ 


(8.28) 

(8.29) 

(8.30) 


and Q d , 
estimate 


(8.31) 
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We may iteratively apply dt to (18.241) to estimate higher-order temporal derivatives; this leads 
us to the estimate 


Qe[q e ] < (1 + PQKeM + ftelc}))X e [f\ + P(* e [v] + X e [c\) \\q e \\l N 
for some universal positive polynomial P with P(0) = 0. Hence 

sup £[q e (t)\< sup (l +P(9\ e [v(t)]+<R e [c(t)]))9\ e [f(t)]) 

0 <t<T 0 <t<T 


(8.32) 


+ sup (1 + P(<R e [v(t)] + 94 e [c(t)])) 
V 0 <t<T 


x exp ( CT 


■f 


V^WaN+ M t )WAN ) dt ) ( \\Qo\\1 N + T 


7 t ii/»ii 

Jo 


AN 


dt . (8.33) 


Similarly, we may derive the bound 

Ozfe] < (l+P(*^[i;]+9^[c]))9^[/|+P(‘H e [t;]+l« c [c]) (||? e ||^ + (m d [v] + ^ d [c])£[q £ ]) ■ (8.34) 
This leads us to the estimate 

[ V[q E (t)]dt < sup (1 + P(9K e [v(t)]+9\ e [c(t)])) [ 9K d [f(t)]dt 
Jo 0 <t<T Jo 

+ sup (1 + P(9t e [u(t)]+SR e [c(t)])) / \\q E {t)\\l N dt 
0 <t<T Jo 

+ [ sup P(D\ e [v(t)] +94 e [c(t)]) J ( sup £[g e (t)]] f (D\ d [v(t)]+D\ d [c(t)]) dt. (8.35) 
\ 0 <t<T J \ 0 <t<T J Jo 

Combining with (18.301) and (|8.33p then leads to the bound 

[ V[q e (t)]dt< sup (1 + P(TK e [v{t)\ + 94 e [c(t)])) f <R d [f(t)\dt 
Jo 0 <t<T Jo 


0 <t<T 


+ sup (P(iHe[u(t)] + 94 e [c(t)])) ( / (*R d [v(t)]+D\ d [c(t)])dt\ sup (Ti e [f (t)]dt) 


+ 


'o 


T sup (1 + P(*R e [v(t)] +iHe[c(t)])) 
0 <t<T 


0 <t<T 


+ sup (P(fR e [v(t)] + £H e [c(t)])) [ {D\ d [v(t)] + D\ d [c(t)])dt 
0 <t<T Jo 


x exp ( CT 


f 


AN 


+ C 


AN 


dt ) ( IkollLv + t 


fwm 


AN 


dt . (8.36) 


8.3.3. Passing to the limit. The estimates (|8.33l) and ()8.36l) provide us with estimates of 

sup £{q e {t)\ + [ V[q e (t)\dt (8.37) 

o <t<T Jo 

that are independent of e. We may then extract a subsequence of e values such that q e converges 
to some q, which by virtue of lower semicontinuity, must also obey the estimates (18.331) and 
(18.361) . Multiplying (18.241) by <p G C} (T x [0,T)) and integrating by parts leads to the identity 

[ [ ~Qs [dttp + K? &w({Kl e v)y) - {Kl e c)y) = f f ipK*/f + f q o ip(-,0). (8.38) 

Jo Jr Jo Jr Jr 

We may then send e —>• 0 to deduce that q is weak solution to (18.191) in the sense of (18.201) . This 
proves the following proposition. 
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Proposition 8.5. There exists a weak solution q to (18.191) obeying the estimates (|8.33l) and 

8.4. Estimates for the solution to (18.11) . We are now ready to return to the problem (18.11) . 

We begin with estimates of the v and c terms that arise in ([8.II) . They are simple variants of 
previous nonlinear estimates (for example Lemmas 13.31 and 16.21 Theorem 15.11 and Proposition 
EH, so we omit the proof. 

Proposition 8.6. Let v and c be given on Ll as in (ED and let v, c denote their bounded 
extensions from T = to T. Then we have the following estimates, where P is a universal 
positive polynomial with P(0) = 0: 

Ke[v\ + ^e[c] < P(£ 0 [r]}) (f[«] + S°[v]) , (8-39) 

Kd[v] + K d [c\ < P{£°[p}) (V[u\ + £°[rj\ + P°[r/]) , (8.40) 

and 

MIn + Man < + £ 0 [v}) + V°[v] + IMlLv+ 1 / 2 ) • ( 8 - 41 ) 

Remark 8.7. The term £°[p] is added onto the right side of (18.401) only because T> a [rj] provides 
no control of rj itself. The rj term in £° [rj\ is enough to make up for this deficit. The term 
IMI 4 JV+ 1/2 added on the right side of (18.41 j) for a similar reason, but in this case the regularity 
demands require more than £° [p]. 

Now we prove that the solution produced by the method of characteristics in (18.121) obeys 
various useful estimates. 


Theorem 8.8. Suppose that u and r / are given and satisfy 

r T 


sup (£[u(t)}+£°[p(t)])+f (/D[u{t)\ + V°[p(t)]\ dt 
0 <t<T v 7 J 0 v 7 


< OO 


(8.42) 


and 


dtp = u ■ Af on E, [nj = 0 on E_, and U- = 0 on E&. (8.43) 

Let q be given by (18.121) . Then q is the unique solution to (18.11) . Moreover, the solution obeys 
the following estimates for some universal positive polynomial P with P( 0) = 0: 


sup £[q(t)]< sup 1 + P(£[u(t)}+£ 
0 <t<T 0 <t<T 


+ sup ( 1 + P(£[u(t)\ + £ 

\ 0 <t<T 


sup (SK e [/(i)]) 
0 <t<T 


^(T)[\\qo\\i N + Tj o WfmiNdt), (8.44) 


and 


[ V[q(t)\dt< sup (l + P(£[u(t)] + £°[p(t)]j) [ <R d [f(t)]dt 
JO 0 <t<T v 7 Jo 

sup (p(£[u(t)} + £°[p(t)})) ( I (v[u{t)]+D°[p(t)}+£ 

0 <t<T v ' \.ln \ 


dt ) sup (V\ e [f(t)]dt) 

0 <t<T 


+ E(T)[\\q 0 \\i N + T \\fmi N dt 


T sup (1 + P(£[u(t)\ + c 

0 <t<T 


+ 


0 <t<T 

where we have written 


sup t (p(£[u(t)] + £ 0 [p(.t)])^j (v[u{t)\ + D°[p(t)] 


+ £ l 


dt 


, (8.45) 


"(T) := exp ( CT J^ P(£[u(t)\ + £°[p(t)]) (^N*)] + V°[p(t)\ + \\p(t)\\l N+1/2 ^j dt ) . (8.46) 
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Proof. The equations in (18.43[) guarantee that the computations in (18.31) and (18.41) are valid, and 
so v satisfies (18. 5|> . First consider T = ff + . Proposition 18.51 yields a weak solution q to (18.191) 
on T obeying the estimates (18.331) and (18.361) on all of T. We call q the restriction of q to T; 
Lemma 18.41 then guarantees that q is the unique weak solution to (18.91) on T, but Proposition 
13.11 guarantees that q coincides with the solution produced by characteristics in (18.121) . The 
estimates (18.441) and (18.451) follow easily from (18.331) and (18.361) and Proposition 18.61 A similar 
argument works for T = fl_. □ 

8.5. Some more useful estimates for (18.91) . To conclude our analysis of the transport prob¬ 
lem we record another a priori estimate for solutions to (18.91) that will be useful in the next 
section. 


Proposition 8.9. Suppose that q is a solution to (E2D satisfying 

sup £[q(t)\ + f V[q(t)]dt < oo. 


Assume also that 


o <t<T 


7= sup (||c(t)|| ofc + |Hf)|| cfc ) < oo 
0 <t<T 


for some 1 < k < AN. Then there exists a universal constant C > 0 such that 


ll9(i)IU<e Cl *ll®IU 

for t € [0,T]. In particular, if qo = 0 then 


oClit-s) 


11/00 


sup \\q(t)\\ 2 k <Te 2C " /T ||/00l|jfcds. 

0 <t<T JO 

Proof. Let a € N 3 with |a| < k. Applying d a to (18.91) leads to the equation 

(dtd a q + v ■ Vd a q + cd a q = d a f — /“, 

1 d a q(t = 0) = d a q 0 


where 


f a =Y^ C a,p (d a ~ P v ■ Vd fj q + d a ~PcdPq) . 


(8.47) 

(8.48) 

(8.49) 

(8.50) 

(8.51) 

(8.52) 


f) <a 


Because of the condition (|8.5I) we may multiply by d a q and integrate to deduce the standard 
energy identity 

d \\d a q\\ 2 0 


= [ (divu-2c)^+ / ( d a f-r)d a q. 
Jr z Jr 


dt 2 

From this and the structure of f a we may easily deduce the estimate 


(8.53) 


d 


* l|2 

«llo 


dt 2 - Mruc'“ 2 \\ da< t\\o + C (\\ c \\c k + IMIc fc ) h\\k \\ da Q\\o +j^ da fd a q (8-54) 

for some universal C > 0. Summing this inequality over all |or| < k, we find (since k > 1) that 


1 „ 

l c llc° + 2 ItIIc 1 


, II 112 

d Nl k 

dt 2 


<Cl\\q\\ 2 k + \\f\\ k \\q\\ k , 


(8.55) 


where C > 0 is another universal constant. The differential inequality (|8.55l) and a standard 
Gronwall argument then imply (18.491) . The bound (18.501) follows from (18.491) and the Cauchy- 
Schwarz inequality: 


jjV' 1 -" < e 2Clt (jf‘ II/WIIJ*) 


< te 2Clt f \\f(s)\\ 2 k ds < Te 2C ^ T f \\f{s)\\ 2 k ds. (8.56) 

Jo Jo 


□ 
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9. Local well-posedness of (11.331) with a± > 0 

9.1. Data: construction and estimates. Our goal now is to deal with the initial data. We 
assume initially that 

n 0 € H 4N , 77o € H m+ 1/2 , v^V*r/ 0 G H m , and q 0 G H 4N . (9.1) 

We may then construct the data for the higher temporal derivatives as in Appendix [Bj This 
process leads us to the following estimate. 


Proposition 9.1. Suppose that £[ 7 / 0 ] < Si/2, where £[%] is given by (1A.10D and <5i > 0 is from 
Theorem \6.5l Then 

2 A 


£ a?«(o) 

3 =1 


2 

4N—2j 


dtq( 0 ) 


4A-2/+1 


+ 


dtV{ 0 ) 


< 


4A-2/+3/2 


Consequently, 


INILa + 11 ^7o 11 4 jv+1/2 p Ikollijv • 

(9.2) 


%<%<%. (9.3) 

Proof. The proof is essentially the same as that of Proposition 5.3 of m and is thus omitted. □ 


9.2. Approximate solutions. Suppose that n, q, rj are given. We then define the forcing terms 
F 1 on and F± on £+ according to 

F 1 [u,q,p\ = -(p+q+d 3 p9)(-Kd t 9d 3 u+u-V A u)~ pV A (ti (p)q)-V A 1l-g{q+d 3 p9)V A 9, (9.4) 

F l[q,v\ = -P' + (p+)q+N + + gr]+M + - 1Z+M+ - (7+ div* ((1 + |V*r? + | 2 ) _1/2 - l)V*r/+) J\f+, 

(9.5) 


and 

~F-[q,v] = - l p '(p)qj N- + M gri-J\f--l'R.}Af-+a- div* ((1 + |V*r/_| 2 ) _1/2 - l)V*r/_) A L. 

(9.6) 

Here the function 1Z is determined by (11.301) . We also define the forcing term / on H via 

f[u, rf\ = - di v A (pu) + Kd t 9d 3 p9 - di v A (d 3 p9u). (9.7) 

Finally we define the density function 

P[q,v] = p + q + d 3 p9. (9.8) 

We define the following functionals for use in the forcing estimates: 

2A-1 .. 2A-1 

[?] := X] | 3s« 


2JooM := 

3 =0 


d{u 


4 AT—27—1 


7=o 


4 AT—27—1 


2 A 


?)ooM : = II ^711 4AT— 1 /2 + Y W^tV 

3 = 1 


4A—2/+1/2 


(9.9) 


and 


2 A 

2)2M := X] 

7=o 


\%u 


2 

4A—2/ 


2 A 

4)M := £ 

7=0 


2 

4 A—27 ’ 


2)2 [v\ ■= \\v\\l N 


2A+1 


+ Y l&tV 


3 =1 


2 

4A—27+3/2 ' 


(9.10) 


For the sake of brevity we will write multiple arguments within brackets to indicate sums; for 
example, %)oo[u,q,rj\ = 2)ooM + 2)oo[<?] + 2)oo [v\- 

Our next two results contain the estimates of the forcing terms used in the problems m 
and 18.11 The proofs are again standard nonlinear estimates and are thus omitted. 
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Proposition 9.2. Let F 1 and F 2 be determined by {u,q,rj) as in (J9.4I) (19.61) . Let JT and Foo 
be given by this choice of F l ,F 2 as in (13.61) and (UEZD- Then we have the estimates 

Foe < ?)oo [q, V\{ 1 + PODoolv])) + P(Voo[u, q, ri])^oo[u, q, rj\ + a 2 || V*t/||| jv _ 1/2 P(2 )oo [v]) (9.11) 

and 

Fi % ?)2 [q, v\( 1 + P{ ?)oo W)) + P{Z)oo[u, q, rj\)fQ 2 [«, Q, v\ + a 2 \\^*v\\l N +i /2 p (2)ooW) (9-12) 
for some universal positive polynomial P such that P( 0) = 0. 

Proposition 9.3. Let f be determined by ( u,rj ) as in (19.71) . Let D\ e be defined by (18. 6R . Then 
we have the estimates 


S»e[/] <?)ocM(l + J P(2)oo[r/]))+P(2)ooW)2)oo[r/] (9.13) 

and 

Kdlf] + 11/llLv < ®M (1 + ^(£°M)) + ^(f°W) (m 4IV+1/2 + V°[v}) (9.14) 

for some universal positive polynomial P such that P( 0) = 0. 

The following proposition allows us to estimate %)oo[ u , q, q] in terms of £ and T> with the 
benefit of introducing time factors. 


Proposition 9.4. We have the following estimates: 


and 


sup 2)oo[«(*)] < £q +T f V[u(t)\dt , 

o <t<T Jo 


(9.15) 


sup 2)oo[g(*)] < So +T 2 sup £[q(t)}, (9.16) 

0 <t<T 0 <t<T 


sup 2Joo[i?(l)] < £o +T I V°[q(t)]dt. 
0 <t<T Jo 


(9.17) 


Proof. The estimates follow easily from the fundamental theorem of calculus and the Cauchy- 
Schwarz inequality. □ 


Now we present the construction of a sequence of approximate solutions. 


Theorem 9.5. Suppose that (uo,qo,Vo) satisfy the compatibility conditions ()B.4I) as well as the 
bound m- Further assume that 

< y, (9.18) 

where T[??o] is given by (1 A. 101) and d'i is given by Theorem, 1 6 . ,71 Then there exists aT± = T±{£q) 
such that if 0 < T < T 4 then there exists a sequence {{u n , q n , q n )}f =0 defined on the temporal 
interval [0,T] satisfying the following three properties. First, ( u n ,q n ,r] n ) achieve the initial data 
at t = 0. Second, for n > 1 we have that 


and 


'p n ~ l d t u n - div^n §_ 4 nu n = F 1 [u n ~ 1 ,q n ~ 1 ,r] n ~ 1 ] in Ll 
d t r] n = u n ■ Af n on T, 

-§ A nU n M n = -a+A^ n AT n + Fl[q n ~ 1,^-1] on £+ 

< -l§ A nu n }J^ n = a-A*rj n Ar n - Fflq 11 - 1 ,^- 1 ] on Z _ 

[u n ] =0 on £_ 

vlf = 0 on Sh 

u n {-, 0 ) = u 0 ,q n {-, 0 ) = 770 , 


d t q n - K n d t e n d 3 q n + div^n ( q n u n ) = f[u 11 , if] in Q X (0, T) 
q n (-,0) = q 0 in LI, 


(9.19) 


(9.20) 
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where F l , F 2 , and f are defined by (19.41) (19.71) and p n 1 = p[q n 1 ,rj n x ] is given by (19.81) . 
Third, we have the estimates 


sup (e\u n {t)\+£ a [rf 
0 <t<T V 


+ / ( V[u"(i)\ + ||p"- 1 J"d“+ 1 u”(<)|| bifi) , + D’fo“(<)]) dt < P,(£g), 

sup ||l)"(t)|| 4 N+1/2 < 

0 <t<T 


sup £[q n {t)} + / V[q n {t)]dt < Pfi£fi), 

0 <t<T JO 

Z[h n ](T) < 6 1 and ±p* < p n = p[q n , V n ] < \p * 
for all n > 1, where Pi for i = 1,2,3 is a universal positive polynomial with Pfi 0) = 0. 


(9.21) 

(9.22) 

(9.23) 

(9.24) 


Proof. We divide the proof into three steps. 

Step 1 - Seeding the sequence 

To begin, we extend the initial data to a triple that belongs to the function spaces necessary 
for the construction of solutions. We combine the data estimates of Proposition 19.11 with the 
extension results of Propositions lF.il IF.21 and lF.4l in order to produce a triple (u°, q°,r] 0 ) defined 
on the temporal interval [0,oo), achieving the initial data, and satisfying the estimates 

POO 

sup£[u°(t)] + / V[u°(t)\dt < P 0 (SS), (9.25) 

t> o Jo 


POO 

sup £[q°(t)\+ V[q°(t)\dt< P 0 (£ q), (9.26) 

t> o Jo 

and 

~ 9 f°° - 

sup^ CT [? ? 0 (t)]+sup||r ? 0 (t)|| , + / V a [r] 0 (t)}dt < P 0 {£g) (9.27) 

t> 0 t> o ' Jo 

for some universal polynomial Pq > 0 with Pq( 0) = 0. 

Step 2 - The iteration procedure 

We claim that there exist universal positive polynomials Pi for i = 1, 2,3 such that Pfi 0) = 0, 
and a > 0 (depending on £q) with the following two properties. First, 


mm{Pi(z),P 2 (z),P 3 (z)} > Pq{z) for all z > 0. 


(9.28) 


That is, each of the coefficients of Pfiz), i = 1,2,3, is bounded below by the corresponding 
coefficient of Po(Z). Second, if T < min{a, T 3 } (where T 3 = T 3 (£{f) > 0 is given by Theorem 
m, and the triple ( u n 1 , q n 1 ,rj n x ) is given, achieves the initial data, and obeys the estimates 

,.T 

sup (£[u n ~\t)] + + / (v[u n -\t)\ + V a [ri n -\t)\) dt < Pfi£S), (9.29) 

0 <t<T ' J Jo ^ ' 

sup ||i? n_ 1 (t)||J iV+1/2 < P 2 (£q), (9-30) 

0 <t<T 1 

and 

sup T[g n_ 1 (t)] < P 3 (£fi), (9.31) 

0 <t<T 

then there exists a triple ( u n ,q n ,r ] n ) that solves (19.191) and (19.201) . achieves the initial data, and 
obeys the estimates 


sup (£[u n (t)]+£ a [r, r 
0 <t<T 


+ £ (v[u n (t) ] + ||p- 1 J^ 2JV+1 ^(t)||^ i)t +P CT [7? n (t)j) dt < P\ (£q ), (9.32) 
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sup ||r? n (t)||4 JV+1/2 < P 2 (fo)> ( 9 - 33 ) 

0 <t<T 

an d 

sup £[q n {t)\ + f V{q n (t)\dt < P 3 (f %). (9.34) 

o <t<T Jo 

The proof of the claim is very similar to Step 2 of Theorem 6.1 in HD. so we will only provide 
a sketch of the idea. It suffices to show that if (19.291) (19.311) hold for some choice of P t , i = 1,2, 3, 
satisfying (19.281) . then ( u n 1 q n 1 rf l ) can be constructed and must satisfy (19.321) (19.341) so long as 
the constants and degree of the Pi are sufficiently large and T < a for some small a. 

The first step is to use Theorem 17.21 to produce a ( u n ,r ] n ) solving (19.191) . For this we must 
verify that the hypotheses *)3(<5i) are satisfied. The hypotheses (17.21) and (17.31) are satisfied by 
assumption. The hypothesis (17.41) follows by combining the estimates of Propositions 19.1119.21 
and 19.41 with the bounds (19.291) (19.311) . The hypothesis (17. 5|) requires that p n ~ l = p[g n_1 , i/ 1-1 ] 
satisfies (13.21) and (15.21) . The condition (15.21) follows trivially from (19.291) and (19.301) . To verify 
CL2D we first estimate 


sup || p n 1 (t) -po\\ LOO < [ || d t p n 1 (t)\\ L00 dt< [ {\\d t q n 1 {t)\\ LOO + \\d t rf l 1 (t)\\ LOO )dt 

0 <t<T Jo Jo 

< T sup ^£[q-~\t)]+£^-\t)} < T JPi(£q) + P 3 {£q) < aJP 1 (£g) + P 3 {£%). (9.35) 

0 <t<T v v 

Then we find that 

sup || p n ~ l {t) - po\\ Loo < ^ (9.36) 

0 <t<T * 

if a is chosen sufficiently small with respect to the Pi, P 3 , £q and p*. Hence (13.21) is satisfied 
by the assumption on po in (12.11) . 

We may thus apply Theorem 17.21 to produce the solution pair ( u n ,rj n ) solving (19.191) . To 
derive the estimate (19.321) we sum (17.151) and (17.171) and again employ Propositions 19.1119.21 and 
IQ and (19.291) (19.311) to estimate the forcing terms. The actual derivation of (19.321) is tedious 
and will be omitted, but we will point out the key observation. The estimates of Proposition 
19.41 guarantee that any appearance of Pi in the resulting estimates is multiplied by at least one 
factor of T and so by choosing a small enough (in particular a bound like a < «o(l + £o) _m f° r 
ao small and m large is needed to reduce the degrees of various polynomials appearing in the 
estimates) and the constants and degrees of Pi large enough, we can show that (19.321) holds. 
The estimate (I9.33|) follows from (17.181) via a similar argument. 

The second step is to use the newly-constructed pair ( u n ,rj n ) to construct q n , the solution 
to (19.201) . through an application of Theorem 18.81 The hypotheses of the theorem are satisfied 
due to (19.321) . (19.321) . and (I9.19|) . The estimates (|8.44|) and (|8.45l) then lead to the estimate 
(I9.34p by employing Propositions 19.31 and 19.41 and arguing as above, except that we use (19.321) 
and (19.331) since the forcing terms are generated by ( u n ,r ] n ). 

Step 3 - Constructing the sequence 

To conclude the proof we combine the previous two steps as follows. We set ( u°,q°,rf) to 
be the triple constructed in Step 1. The bounds (19.251) (I9.27|) imply (19.291) - (19.311) with n = 1 
due to (19.281) . We then set T 3 = min{P 2 ,«} and use Step 2 to construct (u 1 ,^ 1 ,?? 1 ) satisfying 
(19.321) (19.341) and solving (19.191) and (19.201) . The bounds allow us to iteratively apply Step 2 to 
produce ( u n 1 q n ,rf l ) for n > 2. This produces the sequence {(u n , q n , r] n )}'^ = 0 satisfying (19.19D - 
(19.231) for n > 1. It remains only to prove (19.241) . The estimates of p n = p[q n , p n ] can be derived 
exactly as in (19.351) (19.361) . The estimate of £[rj n ](T) can be derived similarly: 


2[v n m<h[r, 0 ] + 3T [ T \\d tri n (t)\\l N _ 1/2 dt<^ + 3T 2 sup £°[ V n (t)} 

1 Jo 4 0 <t<T 

< Ml + 3T 2 Pi(<? 0 ct ) < Mi + 3a 2 P\ (£q) < Si (9.37) 

if a is further restricted. □ 
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9.3. Contraction. We wish to ultimately show that the sequence {(u n , q n , p n )}ff = Q contracts 
in some lower-order regularity space than that given by (I9.21l) - (|9.23l) . Our goal now is to 
prove such a contraction result. We will prove the result in a somewhat more general context 
than within the sequence {(u n , q n , p n )}^=o in order for the result to be applicable in proving 
uniqueness of solutions to (11.331) . 

Before stating the result we define the low-regularity norms in which contraction occurs. We 
define 

2UooM = IM| 2 + II^Ho and W 2 [u] = ||ii||§ + \\d t u\\l , (9.38) 

= ll ? ?ll 5/2 + ll^ r /ll 3/2 + cr ll v * T ?ll 2 + 0 ‘ll v * 5 t 7 ?llo and2U 2 [r/] = cr 2 IMI 7/2 + ||^||?/ 2 > ( 9 - 39 ) 

and 

W 00 [ q } = \\q\\l + \\d t q\\ 2 1 . (9.40) 

Now we state our contraction result. 


Theorem 9.6. Suppose that the triples (u l ,q l ,rf) and (u\p\0) fori = 1,2 satisfy 

'p l dtu l — div^i §_4 iu l = C*] in D 

dprf = u l ■ M l on £ 

-S A iU i M i = -d+A^W' + F 2 \p\ C] on £+ 

- [S^u 4 ] M* = a-A^r ] i Af i - Fl\p\C] on £_ (9.41) 

[td] =0 on £_ 

u l _ = 0 on Sb 

y{-, 0) = u o ,p l {-,0) = r/o, 

and 

(dtq 1 — K l dt6 l d^,q r + di -v A i{q l u l ) = f[u l , rf) in x (0, T) 

\g l (-,o) = g 0 in 

where A\ Af l , 9 l , and K l are given by rf', and F 1 , F±, and f are defined by ()9.41) - (19.71) and 
P % = p[p l ,C] is given by (19.81) . Further suppose that 


(9.42) 


and that 


max < sup sup £[u l ,q t ,p t ], sup sup £ l ] > < M 2 , 

i 0<t<T i 0 <t<T 


£W](T) < Si and < p 1 = pip 1 A 1 } < | P*, 


(9.43) 


(9.44) 


where £ is given by (E3H and (5i is given by Theorem 1 6. FA 

There exist universal constants 7 > 0 and a constant T 5 = Tg(M) £ (0,1) such that if 
0 < T < T 5 and 


max 


sup sup ||r/*(t)|L,sup sup ||C*(i)|L > < 7 2 , 

i 0 <t<T i 0 <t<T 


(9.45) 


then 


sup SBoo [u 1 (t)-u 2 (t),q 1 (t)-q 2 (t),rj 1 {t)-r] 2 (t)]+ / ^[u 1 (t) - u 2 (t), p 1 (t) - p 2 (t)]dt 

0 <t<T J 0 

sup TB 00 lv 1 {t)-v 2 {t),p 1 (t)-p 2 {t)A 1 {t)-C 2 (t)] + \ [ %B 2 [v 1 (t)—v 2 (t),( 1 (t)—( 2 {'t)}dt. 

z 0<t<T * JO 

(9.46) 


Proof. We divide the proof into several steps. 

Step 1 - Differences 

To begin we define u = u 1 — u 2 , q = q 1 — q 2 , p = p 1 — p 2 , v = v 1 — v 2 , p = p 1 — p 2 , and 
f = I) 1 — (f 2 . Then we subtract the equations (19.411) with i = 2 from the same equations with 
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i = 1 to deduce equations for (■ u,r ]). We then apply d a for a € N 1+2 with |a| < 2. This results 
in the equations 


' p 1 dtd a u — div^i §_ 4 i d a u = div^i §>d<*( A i- A 2 )U 2 + d a H l + H 1,a 
dtd a r) + u 2 ■ \7*d a r] = d a u ■ Af 1 + H ’ > 0 

-S^uAf 1 = -a+A^pAf 1 + + d a H\ + H 

< - [S^u] Af 1 = a.A^rjAf 1 + [S^i-^n 2 ] Af 1 - d a H 2 _ - H 2 _! a 
= 0 
d a u _ = 0 

^ Q n(-,0) = 0, d a p(-, 0) = 0. 


in fi 
on £ 
on £ + 
on £_ 
on £„ 
on £;, 


(9.47) 


A similar argument with (I9.42|) but not employing derivatives yields an equation for q: 


dtq — K 1 dtd 1 dsq + div_ 4 i (git 1 ) = H 4 in 12 x (0, T) 
q(-, 0) = 0 in 12. 


(9.48) 


Here we have written the forcing terms as follows: 

H 1 :=F l [v\p\ C 1 ] - F\v 2 ,p 2 ,<?) - (p + d 3 p0[(])d t u 2 + G 1 
Hl-.= Fl\p\?}-Fl\p\e] + Gl 
H 4 ^/[nVj-ZIuVj + G 4 , 

where 0[£] is determined by (, 

G 1 : = div^i__ 42 ) S^u 2 

G\ := —cr + A*r? 2 (A7 1 — AT 2 ) + E A 2U 2 (J\f 1 — Af 2 ) 

G 2 _ := —a-A*r) 2 (Af 1 - Af 2 ) + [S^u 2 ] (Af 1 - Af 2 ) 

G 4 := ((K 1 - K 2 )d t 9 l + K 2 (d t Q l - d t 9 2 ))d 3 q 2 - div^i.^ )(q 2 u 1 ) - div A 2 (q 2 u). 


(9.49) 


(9.50) 


H 1,a ■■= (d a (div A i S A iu) - div^i S A id a u) + (5"(div^i S^i-^u 2 ) - div^i S d a (A i_ A 2 } u 2 ) 
H 2 f := (d a (S A1 uAf 1 ) - S^uAf 1 ) + {d a (S {A i_ A 2 ) u 2 Af 1 ) - S da{A x_ A 2 )U 2 Af 1 ) 

- a + {d a [A*r)Af 1 ] - A^rjAf 1 ) 

H 2 J a := {-d a [S^rnAA 1 ] + [§^i<9“n]] Af 1 ) + (-«9“ [S^i.^ju 2 ^ 1 ] + [S^.^u 2 ] Af 1 ) 

- (t — (<9 Q [A*?7)A/’ 1 ] - A it d a pAf 1 ) 

H 3 ’ a ■= (-d a (u 2 • V*? 7 ) + u 2 ■ V,d a r t ) + (d a (u • Af 1 ) - d a u ■ Af 1 ) . 

(9.51) 


We have written the forcing terms in this manner in order to isolate those terms depending on 
(■ v l ,p\C) from those depending on (u\q\p' 1 ) and to single out some special delicate terms. 
Step 2 - Energy estimate 
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The starting point for the contraction analysis is a basic energy estimate for (19.471) . Arguing 
as in Lemma 16.11 leads us to the equality 


d_ 

dt 


lT i\d a u\ 2 , f \d a V \ 2 , |5“V*p| 2 ' 

P J —n -1- / —^-1- cr- 

4E 


+ 


p J 1 


3Q„.|2 


in 


^id a u\ + J l p |div_ 4 i d a u 


= [ J 1 d a u-{d a H 1 +H 1 ’ a )~ [ d a u ■ {d a H 2 + H 2 ’ a ) 

Jn J e 

+ [ (-5“r/ + (jA,9 Q 77)(u 2 -V*a Q 7 ? -iL 3 ’ a )+ [ d a r](d a u ■ Af 1 ) 

J E J E 


+ 


[ &t{J l p 1 ) 

Jn 


iu\d a u\ 


j f -7r°A ldau : ^°d°‘{A 1 -A 2 ) u2 + ^V(div 4 i d a u)( div a c«(^i__ 42 ) u 2 ). 

(9.52) 


Let us define 


m= £ ghill! +g ia°v.^)i 2 ' 


aeN 1+2 

|a|<2 


and 


®(t)= j; na“«(*)||f. 


(9.53) 


(9.54) 


aeN 1+2 

| a |<2 


We sum (19.521) over a € N 1+2 with |a| < 2; applying Proposition IE.31 and arguing as per 
usual (as in Lemmas 13.31 and 16.21 Theorem EH and Propositions 16.41 and 19.21) to estimate the 
various nonlinearities, we derive the differential inequality 


j t m+cn(t) < p( 7 )v^R\/^ icl 

+ (1 + P(M))V2HooM (V2»oob,P,C] + V2Boo[^,C]) + (1 + P(M))^W 2 [uWW 00 \p 1 C\ 

+ (1 + P(M))2U 00 [u, rj\ + (1 + P{M))y/ 2Uoc [u, rj\ \J 2U 2 [u, r/\ (9.55) 

for some universal positive polynomial with P(0) = 0 and a universal constant C > 0. We 
should note that two of the terms appearing on the right side of (19.521) require some delicate 
treatment. The first are terms involving Vp in H 1 . In order to handle these when two horizontal 
spatial derivatives are applied, we must integrate by parts to move one horizontal derivative 
onto J l d a u and reduce to only two derivatives on p, which is all that is controlled by SU^p], 
The second are terms involving a± multiplying two spatial derivatives of i n P±; these give 
rise to the term P(y)y^2U2 [u] ^2X12 [C] • The key part of this is P(y), which appears because the 
nonlinear terms only involve spatial derivatives of Q. 

Integrating (I9.55P in time, using the fact that 11(0) = 0, and applying the Cauchy-Schwarz 
inequality then yields the bound 


sup it(f) + [ 93 (t)dt<P(j) f (2U 2 [C ]+W 2 [u])dt 
<t<T Jo Jo 

+ Vt(1 + P(M)) [ (2H 2 [u,C]+2 V 2 [u,r t ])dt 
Jo 


+ (VT + T)(l + P{M)) sup 0mo>,p,C]+2IJooM]), (9.56) 

0 <t<T 

where again P is a universal positive polynomial such that P(0) = 0. 

Throughout the rest of the proof we will let Z denote a quantity of the form 


Z ~ (1 + P(7)) (RHS of dHESD) , 


(9.57) 
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where P is some universal positive polynomial such that P(0) = 0. From one estimate to 
another the polynomials and constants may change, but the structure of Z does not. 

Step 3 - Improved u estimates 

The usual trace theory allows us to estimate 

/ IMIhV 2( E ) ~ H <9 "' u lli~ [ W(t)dt<Z. (9.58) 

J ° J ° J ° 

\ a \<2 

We may apply the elliptic estimate of Proposition IE. 41 which is applicable due to the first 
estimate in (19.441) . to bound 


W\l ~ l|p 1<9 HI^ + lldiv^i S(^1-^2)U 2 + H 1 ^ + ||«||^5/a (E) 


All 2 


~ \\P llL° 


? + ||Vp 1 ||L 


o + (1 + P(M)) C] + 2UooN) + |M|tf5/2 (S) 


sup 
0 <t<T 


0 <t<T 


< Qj + (1 + P(M))2Hoo [u] + (1 + P(M)) (Wnfap, c] + SUooM) + ||n||^5/2 (s) . (9.59) 

Here in the third inequality we have used the second estimate in (|9.44l) . Hence 

/ \\u(t)\\ldt< [ T V(t)dt + Z<Z. (9.60) 

Jo Jo 

We improve the L°° in time estimate for u by employing Lemma lD.il 

sup ||u(t)||| < /" (||u(t)||| + \\dtu(t)\\l)dt < f (||«(t )||3 + 9J(t)) dt <Z. (9.61) 

0 <t<T Jo Jo V J 

Combining (19.561) . (19.601) . and (19.611) then provides us with the bound 

sup 2B 00 ['u]+ [ W‘ 2 [u]dt < Z. 
o <t<T Jo 

Step 4 - Improved q estimates 

Now we improve the estimates for q. Note first that we already have the bound 
\l + \\9tv(t)\\ 2 0 + o ||V*77(t)||| + cr ||V*9 f r ? (t)|| 2 ) < sup it(t) < 2. 


(9.62) 


(9.63) 


By solving the second and third equations in (I9.47|) for oAq and employing (I9.60|) . (19.611) . and 
(19.631) . we may then bound 


\\v(t)\\y 2 dt < (\\q(t)\\l + \\aA*q(t)\\l /2 ^ dt 

<T sup \\q(t)\\l+ [ ((1 + P(m))W 2 [u] + P('y)W 2 [C]) dt 
0 <t<T Jo 

+ T(1 + P(M)) sup (W 00 lp,(]+W 00 [rj\)<Z. (9.64) 

o <t<T 

Next we employ the kinematic equation in (19.471) along with the transport estimates of Propo¬ 
sition 2.1 of [7] to bound 

Il^(^)|l5 /2 < exp ( CT J T f 0 h'N 1 (.1)111/2 dt - ( 9 - 65 ) 

Hence, if we assume that T 5 M < 1 we may bound the exponential term above by a universal 
constant and then estimate 

s^P lk(i)ll5/ 2 < r ( 1 + P (7)) / \\u(t)\\l dt<TZ. (9.66) 

0 <t<T Jo 

Then we solve for dtq in (19.410 to bound 

sup \\d t q(t)\\l/ 2 < M sup \\q(t)\\l /2 + (1 + P{l)) sup \\u(t)\\ 2 2 < MTZ + Z < Z (9.67) 

0 <t<T 0 <t<T 0 <t<T 
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since MT < MT§ < 1. Similarly, we solve for to estimate 

[ \\ d tV(t)\\ 2 l/2 dt ~ MT SU P 11^(0111/2 + (! + 7) [ \\d t u(t)\\ldt 
Jo ' 0 <t<T Jo 

+ MT sup (||u(t )||2 + ||^( 0 lll/ 2 ) ~ ( 9 - 68 ) 


0 <t<T 

Summing (I9.63p . (19.641) . and (|9.66|) — (|9.68|> then yields the bound 

sup 2Hoo[r/] + [ 2 U 2 [q]dt<Z. 
o <t<T Jo 

Step 5 - Estimates of q 

Next we employ Proposition 18.91 to get estimates for q. First we find that 

sup |k(Olll<exp(C'(l + P(M))T)T f \\H\t)f 2 dt. 

0 <t<T Jo 


(9.69) 


(9.70) 


If we further restrict T§ so that (1 + P(M))T < 1 we can again treat the exponential as a 
universal constant. Then 

sup \\q(t)\\ 2 2 <T [ \\H 4 (t)\\ 2 2 dt<T(l + P(M)) f (2D 2 [«] + r,}) dt 


0 <t<T 


T(1 + P{M))Z. (9.71) 


Next we use (19.411) to solve for dtq and estimate 


sup \\d t q{t)\\l < P{M) sup \\q(t)\\l+ sup ||P 4 ||^ (9.72) 

0 <t<T 0 <t<T 0 <t<T 

The term H l may be estimated as follows. First we use the fact that u(t = 0) = 0 and 
r](t = 0) = 0 to estimate 

j~T 

sup \\dtv(t)\\l /2 + sup ||Ti(t)Hi<r / +\\d t u(t)\\f\dt<TZ. (9.73) 

0 <t<T 0 <t<T J 0 V ' J 

Then we bound 

sup \\H\t)\\ 2 < P(M) sup \\rj(t)\\l /2 +P(M) sup \\d t r](t)\\ 2 1/2 
0 <t<T 0 <t<T 0 <t<T 

+ P(M) sup IK0II? + (1 + P(t)) sup \\u(t)\\l 

0 <t<T 0 <t<T 

<TP{M)Z + {l + P( 1 ))Z<Z (9.74) 


if T 5 is further restricted so that T§P(M) < 1. Note here that we have crucially employed the 
T factor appearing on the right side of (19.661) and (19.731) . We can now combine (19.711) . (19.721) . 
and (19.741) to deduce that 


sup SUoo [q] < T{ 1 + P(M))Z + Z<Z (9.75) 

0 <t<T 


if T 5 is further restricted. 

Step 6 - Synthesis 

Now we sum the estimates (19.621) . (19.691) . and (19.751) to deduce that 

sup 2Boo [u,q,rj](t)+ f 2 V 2 [u,rj}(t)dt < Z. 
0 <t<T Jo 


(9.76) 
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Assuming that 7 and T 5 are sufficiently small and using the previous inequality, we may absorb 
the terms involving (u,q, p) from the right side to the left; this results in the estimate 

sup 2Uoo [u,q,p](t) + [ W 2 [u,q\(t)dt<P(p/) [ 2H 2 [C ]dt 

0 <t<T JO Jo 

+ Vt {1 + P(M)) [ W 2 [v,(\dt + {VT + T)(l + P(M )) sup (9.77) 

Jo 0 <t<T 

By further restricting 7 and T 5 we deduce that (19.46|i holds. □ 


9.4. Local well-posedness. We now have all of the ingredients necessary to prove a more 
general version of Theorem 12.11 in the case a± > 0. 


Theorem 9.7. Assume that a± > 0. Suppose that (uo,qo,Po) satisfy the compatibility condi¬ 
tions (El as well as the bound m, and that 

< f, (9-78) 

where T[?/o] is given by (IA.10I) and <5i is given by Theorem 1 6 . 51 Further assume that 



(9.79) 


where 7 is as given in Theorem 1.9.61 Set Tq = Tq{£q) = min{T 4 (£g ), Tq(Pi{£q ) + -^M^o )}, 
where T 4 and P± , P 3 are given by Theorem 1.9.51 and T 5 is given by Theorem 1.9.51 

If 0 < T < Tq then there exists a triple ( u,q,r /) defined on the temporal interval [0, T] 
satisfying the following three properties. First, ( u , q , 7 ) achieve the initial data at t = 0. Second, 
the triple uniquely solve (11.331) . Third, the triple obey the estimates 


sup 
0 <t<T 


(£[«(()] +£>(()]) 


V (“(*)]+ pJBf N+ 1 u(t) 


sup 11^)11^1/2 <P 2 (^), 

0 <t<T 


dJ<Px (£%), 

(9.80) 

(9.81) 


and 


sup £[q(t)\ + f V[q{t)]dt < P 3 (^o )> ( 9 -82) 

0 <t<T Jo 

£[ri}(T) < 61 and < p = p[q,rj\ < ^p*. (9.83) 


Proof. The proof is very similar to that of Theorem 6.2 in HU, so we will only provide a quick 
sketch. 

First we use Theorem 19.51 to produce a sequence of approximate solutions {(u n , q n , q n )} < ?f=o 
on the temporal interval [0,T]. The uniform estimates (19.21 1) (I9.24j) along with standard com¬ 
pactness and weak compactness arguments yield a subsequence converging to a limiting triple 
(■ u,q,r]) that achieves the initial data and satisfies the estimates (19.8011 (19.8311 . 

Because we only know the convergence of a subsequence, we cannot immediately pass to 
the limit in (19.191) and (19.201) . Instead we first use Theorem 19.61 to deduce that the sequence 
{(u n , q n , rj n )}ff = 0 actually contracts in the lower regularity norm defined by (19.461) . In order 
to apply the theorem we must verify that (19.441) and (19.451) are satisfied; these follow from 
(19.781) and (19.791) . an argument like that used in (19.35H (I9.36H . and a further restriction of time. 
This low regularity convergence, when combined with the bounds (19.211) (I9.24H and various 
interpolation arguments, shows that the original sequence actually converges to (u, q, if) in a 
regularity class slightly larger than that defined by (19.211) (19.241) but more than sufficient for 
passing to the limit in (19.191) (19.201) . We deduce then that (u, q, p) satisfy (11.331) . The uniqueness 
claim follows from another application of Theorem 19.61 □ 
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10. Local well-posedness of (J1.33I) with a± = 0 
We now state a result on the local existence of solutions to (11.3311 without surface tension. 


Theorem 10.1. Assume that a± = 0. Suppose that (uq, go, ?/o) satisfy the compatibility condi¬ 
tions (IB.411 as well as the bounds (12.Ill , and that 


< y, 

where £[??o] is given by (1A.10I1 and <5i is given by Theorem 1 6. .51 Further assume that 



( 10 . 1 ) 


( 10 . 2 ) 


where 7 are as given in Theorem, 1.9.61 Set Tq = Tq(£q) = min{T 4 (£g), T^{P\{£q ) + P^{£q)}, 
where T 4 and P ±, P 3 are given by Theorem 1.9.51 and T 5 is given by Theorem, 1.9.51 

If 0 < T < Tq then there exists a triple (u, q, 77 ) defined on the temporal interval [0, T] 
satisfying the following three properties. First, ( u , q , 7 ) achieve the initial data att = 0. Second, 
the triple uniquely solve (11.3311 . Third, the triple obey the estimates 


sup 

0 < t<T 


(^£[u(t)] + £°[r](t)}^ 


rT / 2 

T>[u(t)\ + pJdf N+1 u(t ) 


(0 HA* 


sup \\vml N+ i/2<M£S), 

0 < t<T 


+ V°[r,(t)})dt<P 1 (£S), 

(10.3) 

(10.4) 


and 


sup £[q{t)\ + f V[q{t)]dt < -P 3 P 0 ), (10.5) 

0 < t<T JO 

£[v)(T) < 5i and i p * < p = p[q,rj\ < ^p*. (10.6) 


Proof. The proof follows from an argument similar to that used in the proof of Theorem 19.71 
but actually somewhat easier and more akin to that used in Section 6 of m- We will provide 
only a sketch of the ideas. 

The main difference between the method to produce solutions to (11.3311 with a± > 0 and the 
method used with <j± = 0 lies in the use of the k approximation, which replaces the kinematic 
transport equation for p with a parabolic problem. This is essential in studying the problem 
with surface tension, as it leads to a regularity gain for 77 that enables us to treat cr±A* 77 ± as a 
forcing term when solving for u. However, when a± = 0, this regularity gain is unnecessary, as 
only 77 appears as a forcing term in the u equation. In place of the parabolic problem we simply 
study the kinematic transport problem directly, using Theorem 5.4 of m to produce solutions 
and derive estimates. 

We then proceed essentially as in Section [9) First we prove that Theorem 19.51 holds with 
a± = 0. The iteration scheme begins with a triple (u”" 1 , q n ~ l , rf 1 ^ 1 ) and then uses Theorem 
13.41 to produce u n . Then the equation dtp 71 = u n ■ M n is solved using Theorem 5.4 of [11]. 
Then (u n ,p n ) are used to solve for q n in (19.201) by way of Theorem 18.81 Next we observe that 
Theorem 19.61 remains true as stated with a± = 0. Finally, we combine these two theorems to 
produce a sequence {(u n , q n , p n }^ =1 of approximate solutions that remain uniformly bounded 
at high regularity and contract in a lower-regularity norm. □ 


Appendix A. Energy and dissipation functionals 


Here we collect the definitions of various functionals that are used throughout the paper. We 
define the energies associated to (u, q, 77 ) via 

2 N 


|p 

3= 0 


u 


2 

41V—2j 


(A.l) 
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2N 


%] = MIn + J2\\ d t q 

3 =1 


4AT-2j+l 


and 


2 N 

£ a iv\ = J2\\ d t r] 

3=0 


4N-2j 


+ cr 


<9/V*r/ 


4Af-2j 


We define the corresponding dissipation functionals via 

2N 

v \ u \ = ^2\\ d i u 

3=0 


and 


27V+1 


V a [i]\ = a 2 \\v\\l N+3 /2 + H^llLv-i + I ^t 1 ! 

3 =2 

For r] we also need to define some improved terms: 

2N 

£ a [v\ = £ a [v\ + J2\\ d t r i 


4N-2j+2 


3 =1 


4Af-2j+3/2 


(A.2) 


(A-3) 


4AT-2j+l 

(A.4) 

27V+! 2 


v 

1 4AT-2j+2 

(A.5) 


(A.6) 


(A.7) 


and 


2N 


V a W\ = V a [iq\ + \\d t ri\\l N _ 1/2 + ||^r? 

J=2 


4IV—2j+5/2 


We must also define the term 


£[t/] (T) = SUp 1127(7) 1147V—1 /2 
0 <t<T 


(A.8) 


(A.9) 


In estimating this term we often refer to the following term associated with the data: 

^[ 770 ] = H 770 II 4 JV— 1/2 ■ (A.10) 

For the data ( 770 ,( 70 , 770 ) we define 

£0 = II u 0 II 4 JV IkolUiv + ll 7 ?© II 4 JV+I /2 + a 11V * 770 114JV > (A.11) 

and when {(c^u(0), djq(0), c^77(0))}|= 0 are known we write 

2N 2N „ 


3=0 


£0 M=Xj3t“(°) ^o[®] = lk(o)||Lv + 5Z N 9 (°) 

z —' N 4N—2j z — 4 II 


3 = 1 


4JV-2.J+1 


2N 


and £0 [ 770 ] = 1177(0)11^4/2 + cr 1^*77(0)11^- + ||^77(0) 

i=i 


4N-2j+3/2 


(A.12) 


We will often abbreviate 

£q = 4>o] + #[*)] + £q [t?o] an d £q[ u 0 j 770 ] = £q[uq] + £q [ 770 ]. 


(A.13) 
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Appendix B. Compatibility conditions 


Here we record the system of compatibility conditions that the initial data ( uo,qo,rjo ) must 
satisfy in order to produce high-regularity solutions to (11.331) . To state the compatibility condi¬ 
tions we must first show how to construct (d{u(-, 0),dj q(-, 0), cP t r](-, 0)) for j = 1,..., 2N from 
the triple (u 0 ,q 0 ,r] 0 ). 

For the purposes of constructing the temporal-derivative data we rewrite the first, second 
and third equations in (11.331) in the form 

d t Ti = Fi (u, 7?) 

dtq = F 2 {u,q : ri,dtri) (B.l) 

dtu = F 3 (u,q,r],d t r]). 


Assuming that we are given {(d^u(-,0),d^q(-,0),d^r](-, 0))}], =0 for some j € {0,... , 2N — 1}, we 
construct (c^ +1 ti(-, 0), dj +1 q(-, 0), 0)) as follows. First we apply d{ to the first equation 

in m and define 


d} +1 v(-, 0 ) = dlFi(u,ri) 


(B.2) 


which is possible because all terms appearing on the right are already known. We can now 
perform a similar operation on the second and third equations in (IB. ID . setting 


d{ +1 u(-,0) = d J t F 2 (u 1 q,T],dtr ]) 


and d{ +1 q(-,0) 


diF 3 (u,q,Ti,d t Ti ) 


(B.3) 


both of which can be computed in terms of known quantities since we have already computed 
di +1 rj(-,0). Using this argument, we may inductively define {(d J t u(-,0),d{q(-,0),d{ri(-,0))}™ 1 
as desired. 

We may now state the compatibility conditions. We say that (uo, qo, Vo) satisfy the compat¬ 
ibility conditions at level 2N if 


'&{ {P\p)qM - Sa(u)AT) | t =o = d{ 
{ p \P)qN ~ §a(uW) |t=o = 
dj/u | t=0 = 0 
dfu- t =o = 0 


(pigrjAf - v+'H+N ~ FFf ) | t=0 
di (IpI + cr-P-FT - pi} A0 


t =0 


on E + 
on E_ 
on E_ 
on Eb 


(B.4) 


for j = 0,..., 2N — 1. 


Appendix C. Poisson extension 


We will now define the appropriate Poisson integrals that allow us to extend q±, defined on 
the surfaces S-t, to functions defined on fi, with “good” boundedness. 

Suppose that E + = T 2 x {£}, where T 2 := ( 2ttLiT) x ( 2nL 2 T ). We define the Poisson integral 
in T 2 x (—oo,£) by 




E 


e ii-x’ 

2lTyJ L\L 2 




(C.l) 


where for £ E (L 1 1 Z) x (L 2 1 Z) we have written 


m 



/(*0 


e ~i^ x ' 

2tx \/ L\L 2 


dx r . 


(C.2) 


Here ” stands for extending downward and stands for extending at x 3 = £, etc. It is 
well-known that : F[ s ( E + ) —> H s+1 / 2 ( T 2 x (—oo, £)) is a bounded linear operator for s > 0. 
Certain improvements of this are available when we restrict to H; we refer to the appendix of 
m for details. 

We extend r/ + to be defined on H by 


r]+(x', x 3 ) = V + r] + (x',x 3 ) := V-/r] + (x', x 3 ), for x 3 < 1. 


(C.3) 
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If 77 + € H s 1 / 2 (S_|_) for s > 0, then g + € H S (Q). 

Similarly, for E_ = T 2 x {0} we define the Poisson integral in T 2 x (—oo,0) by 

V-,o fix) = V —^==e^f{f). 


«e(Lr 


)x(Li 1 Z) 


‘I'K\J L\Li 


(C.4) 


It is clear that V-$ has the same regularity properties as V This allows us to extend t/_ to 
be defined on I2_. However, we do not extend ??_ to the upper domain by the reflection 
since this will result in the discontinuity of the partial derivatives in x 3 of the extension. For 
our purposes, we instead to do the extension through the following. Let 0 < Ao < Ai < • • • < 

A m < 00 for m G N and define the (m + 1) x (m + 1) Vandermonde matrix V(Ao, Ai,..., A m ) 

by V(Ao, Ai,..., X m )ij = (—A j) 7 for i,j = 0,...,m. It is well-known that the Vandermonde 
matrices are invertible, so we are free to let a = (ao, or, •.., a m ) T be the solution to 

V(Ao, Ai, • • • , A m) & — Qm , Qm — (1, 1, ■ • • , I) • (C.5) 

Now we define the specialized Poisson integral in T 2 x (0, oo) by 


V+,of(x) 


E 


27T\/ L1L2 


m 


j =0 


(C.6) 


It is easy to check that, due to f|C.5|> . d l 3 V+fif{x', 0) = d l 3 V-fif{x',0) for all 0 < l < m and 
hence 

d a V+ t of(x', 0) = d a V- t of(x', 0), Va € N 3 with 0 < |a| < m. (C.7) 

These facts allow us to extend rj- to be defined on by 


r]~(x\x 3) = V-r)-(x',x 3 ) 


’P+,oV-(x',x 3 ), x 3 >0 

V- t oV-(x',x 3 ), x 3 <0. 


(C.8) 


It is clear now that if 7?_ G LP~ 1,/2 (E_) for 0 < s < m, then f\- G H s (Li). Since we will only 
work with s lying in a finite interval, we may assume that m is sufficiently large in (IC.5I) for 
fj- G H s (Tl) for all s in the interval. 


Appendix D. Estimates of Sobolev norms 

Here we record an estimate involving space-time norms. 

Lemma D.l. Let T denote either E or fL Suppose that f G L 2 ([0, T]; H Sl (T)) and dt( £ 
L 2 ([0, T]; H S2 (T)) for si > s 2 > 0. Let s = (s 1 + s 2 )/2. Then ( G ^([O, T]; H S (T)) (after 
possibly being redefined on a set of measure 0 ), and 

iiciiI-h, < wcml* + c \\ah H n + c 11&CH W* (d.i) 

for some universal constant C > 0. 

Proof. This is a slight variant of Lemma A.4 of mi that follows from the same argument. □ 

Appendix E. Some estimates involving the geometric terms 

E.l. Coefficient estimates. Here we are concerned with how the size of r] can control the 
“geometric” terms that appear in the equations. 

Lemma E.l. There exists a universal 0 < 5 < 1 so that if H17II5/2 — then 

II^ — l|lL°°(n) 4“ II^IIl°°(D) 4" ll-®lli°°(n) — F’ 

l|4V" — 1 |Il°°(T) 4“ ll-K" _ l|li,°°(r) — am ^ 

11-^ IIl°°(g) 4“ II-^IIl 00 ^) ^5 1- 


(E.l) 









COMPRESSIBLE VISCOUS SURFACE-INTERNAL WAVES 


55 


Also, the map 0 defined hy (11.181) is a diffeomorphism, and 


for all p G L 2 (fl). 



J M 2 <2 [ M 2 

J n 


(E.2) 


Proof. The estimate (IE.II) is guaranteed by Lemma 2.4 of [9]. The estimate (IE.21) then follows 
trivially from (IE. ID . □ 


E.2. Korn’s inequality. Here we record a version of Korn’s inequality that is needed through¬ 
out our analysis. First we record a version involving only the deviatoric part of the symmetric 
gradient, D°, defined by (11.31) . 

Proposition E.2. There exists a universal constant C > 0 so that 

||u||? = IKII? + ||u_||? < C(||D°u + || 2 + ||D°u_|| 2 ) (E.3) 

for all u± G with JuJ = 0 along E and u_ = 0 on E&. 

Proof. The proof is based on the “deviatoric Korn inequality” of Dain, Theorem 1.1 of [6]. For 
details of the proof see the appendix of [H]. □ 

Next we extend Proposition IE.21 to O^. 

Proposition E.3. Assume that p > 0 and p! > 0. There exists a universal 0 < 5 < 1, smaller 
than the 5 appearing in Lemma I E. 1\. such that if 11 11 §/ 2 < then 

IMIi ~ f Q 2 J l B ^ n | 2 + L i ' J l div - 4 u \ 2 ~ Nil ( E - 4 ) 

for all u± G H 1 (Ll±) with [uj = 0 along E and u_ = 0 on E&. 

Proof. Let 5 be as small as in Lemma IE. II With Proposition IE.21 in hand we may then argue 
as in Lemma 2.1 of HH, further restricting 6 as needed, to derive (IE.4D . □ 


E.3. Elliptic estimates. Here we record elliptic estimates for the two-phase geometric Lame 
problem with Dirichlet boundary conditions: 


— div^ §_4U = G 

in H 

u + = h + 

on E + 

u + = U- = h- 

on E_ 

u- = 0 

on Efo. 


Our elliptic regularity result is contained in the following. 


(E.5) 


Proposition E.4. Let k > 4 he an integer and suppose that rj G H k+1 / 2 . There exists 0 < Jo ^ £ 
(where <5 is given by Proposition I E. .H) so that if ||j 7||?_ 1/ / 2 < <5o, then solutions to (IE.5D satisfy 

Nlr< (l|GIU + INU /2 ) (e.6) 

for r = 2,..., k, whenever the right side is finite. 

In the case r = k + 1, solutions to (IE.5D satisfy 

IMU+i ~ (ll G llfc-i + INI/c+1/2) + \\v\\k+i/2 (ll El l | 2 + IHI7/2) ( E - 7 ) 

whenever the right side is finite. 


Proof. Estimates of the form (IE.6D for all r > 2 follow from the standard elliptic regularity 
theory for the problem 


— div §u 

= G 

in H 

u_)_ = h + 


on E + 

u + = u_ 

= h_ 

on E_ 

U— = 0 


on E;,. 


(F.8) 
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With these estimates in hand, we may argue as in Section 3 of [31] to deduce (lE.Gji and (IE. 711 
under the smallness assumption ||??||^_ 1 /2 — ^o- 


□ 


Appendix F. Data extension results 


Next we need some extension results. Our first one allows us to take data with parabolic 
scaling and extend them to space-time functions with the same scaling. The proof can be found 
in Lemma A.5 of mi and is thus omitted. 


Proposition F.l. Suppose that for j = 0, ...,2 N we have that cP t v( 0) € H iN 2j (Q). Then 
there exists an extension u, achieving the initial data, such that 


Moreover, 


sup 

j=Q 


2-A/" 2 /*oo 2iV 

J0 3=0 


AN—2j+l 


for j = 0,. 

.,2N. 

(F.l) 

2N 

2 


3=0 

4N-2j 

(F.2) 


We also need a version of Proposition IF. II that works to extend r/. The difference between r/ 
and u is that Proposition 19.11 provides higher regularity for the time derivatives that we would 
get from strict parabolic scaling. Nevertheless, Lemma A.5 of mi may be readily modified to 
construct an extension. 


Proposition F.2. Suppose that £(0) € iL 4iV+1 / 2 (E) and that y^tTV^O) e H 4N (T,). Further 
suppose that for j = 1,...,2 N we have that cP t Q( 0) € H 4N ~ 2 i +3 ^ 2 (T,). Then there exists an 
extension £, achieving the initial data, such that 

C € L°°([0,oo);iL 47V+1/2 (S)) n L 2 ([0, oo); Pf 4Ar+1 (E)), 

€ L°°([0,oo);iL 47V (S)) nL 2 ([0,oo);77 4Ar+1 / 2 (S)), 
dtC € L°°([0, oo); fL 47V_ly/2 (S)) n L 2 ([0, oo); iL 4Ar (E)), (F.3) 

d{( € L°°([0,oo);iL 47V - 24+3/2 (S)) n L 2 ([0, oo); H m ~ 2j+b / 2 (T)) for j = 2,... ,21V, 

<9 27V+1 C € L 2 ([0,oo);L7 1/2 (S)). 

Moreover, we have the estimates 


2 N 


SU P IICWII 4AT+1/2 +J2 New 


t> 0 


3 = 1 


4AT—2j+3/2 


+ 


2JV+1 

4W+1 + l|S<C(i)llL+ E k'C(t) 

3 =2 


4JV-2j+5/2 
2N 


< 

r^j 


iic(o)iiL +1/2 + E k«o) 

3=1 


dt 


4N-2j+3/2 


(F.4) 


and 


roo II 2 

super || V*C||;Lv+ / cr\\\7*((-,t)\\l N+1/2 dt<a\\\7*((0)\\l N + crY^ dtC(O) • (F.5) 

t> 0 Jo 11 4Y-2J+3/2 


Proof. For each j = 0,, 2 N let tpj € C^°(R) be such that ipf\ 0) = $j,k for k = 0,..., 2 N 
(here ( k ) is the number of derivatives and 5j^ is the Kronecker delta). For the same j let 
fj = ^C(O). Then / 0 € iL 47V+1 / 2 (E), ^V*/ 0 € iJ 4iV (S), and fj € tf 4 iV-2j+3/2( S ) for ^ = 
l,---,2 N. 
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We will construct £ as a sum £ = Y^j=o^j- To construct the Fj we must break to cases, 
considering j = 0 and j > 1 separately. We begin with the latter case, in which case we use the 
Fourier transform given by ()C.2I) for the construction. We define Fj via its Fourier coefficients: 

m,t)=vAt(a) 2 )m (o~ 2j , 


(F.6) 


where (£) = \J 1 + |£| 2 . 

By co 
estimate 


By construction, c£ fc Fj(£, t) = ipj (t (£) 2 )/j(£) (£) 2 ^ fe and hence d^F(-,0) = 6j,kfj■ We 




47V—2fc+3/2 


E«> 


2(47V—2fc+3/2) 


<pf ] (t(o 2 ) 2 m) 2 (o 2{2k - 2j] da 


so 


Similarly, 


f 


< 


EM fc) (^> 
2 


/,(0 (0 


ii/,-ii 47V-2J+3/2 ’ 


2(47V-2j+3/2) 




sup 

f>0 




< 


47V—2fc+3/2 


II /,' 1147V—2j+3/2 for k — 0,..., 2N. 


(F.7) 


(F.8) 


d?Fj(;t) 


dt= \ 

47V—2fc+5/2 ,7 0 


F?\t(o 2 ) 2 m 2 (o 2{2k ~ 2j) dt 


^2(47V-2fc+5/2) 

£ 

> i 9 9 

Ekf’iu?) 2 ) />(0 e) 2( “- 2i+5/2) <it 


o / /*oo 

2(47V-2j+5/2) 


£ 

E|/,(e) 2 ( 0 2(4JV - 2j+5/2) 


(£) 2 

^• fc) (r) 


/o 


(it 


1 ( fc ) ^ 


dr 




(k) 


L 2 (R) 




2(47V—2j+3/2) 


< 


~ 7) 


(*) 


L 2 (R) II^TV—2J+3/2 ’ 


SO 


L 


00 2 


(F.9) 

(F.10) 


m _ 2k+5/2 dt ~ HZ? H 47 V—2J+3/2 for k = 0, . . . , 27V + 1. 

It remains to handle the case j = 0. We define Fo(£, t) = <Po{t (£))/o(£)- Then we may argue 
as above to deduce the bounds 

/ 27v „ \ 


sup \\Fo(-,t)\\l N+1/2 + J2\\ d t F o(-,t) 


t> 0 


k =1 


47V—2fc+3/2 


< 
rv-* 


II /01147V-I-1/2 


and 


roo / 2A^+1 

/ l|To(-,t)||^ +1 + ||9 f T 0 (-,t)||2 Ar + £ 5fF 0 (-,i) 

^ V 77=2 


47V—2fc+5/2 


dt < 


47V+1/2 


for all k > 0. We may also estimate 

12 


sup||V*Fb(-,f)|| 4JV + / || V*Fo(-, t) 114jv+ 1/2 dt < IIV */o || 4 jv ■ 

t>o Jo 


(F. 11) 


(F.12) 


(F.13) 



























































58 


JUHI JANG, IAN TICE, AND YANJIN WANG 


Now, since we set f = Ylj=o we ma y sum (EH), (IF. 101) . (IF. Ill) , and (IF.121) to deduce that 
(1F.4[) holds. To prove (IF.51) we sum (IF.81) and (IF.101) with k = 0 and (IF. 131) : 

POO 

Super ||V*CM)llLv + / CT ||V*CM)llLv+l/2<*f ~ Super ||V*Fi(-,t)||4 Ar 

t >0 Jo t >0 

poo 2 AT /»oo 2N 

+ ^l|V*Fo(-,t)||4 Ar+1/2 dt + suprj £llv,-fX-,i)llb + f a^||V.^(.,<)||^ +I/2 


2N 


< 


G 


iv.e(o)iiiv+<rEF'«o) 4K _ 2 . +3/2 - (F.i4) 


3 = 1 


□ 


We will also need the following simple variant of Proposition IF. 21 We omit the proof for the 
sake of brevity. 

Proposition F.3. Suppose that 5(0) € i? 4Ar_1 (S) and that for j = 1,..., 2 N we have that 
<9;?5(0) € Ff 4Jv ^ 2 - 7_1 / 2 (S). Then there exists an extension 5, achieving the initial data, such 
that 

dl~ € L°°([0, (X)); n L 2 ([0, oo); H iN ~ 2 ^ (£)) for j = 1,..., 2 N. (F.15) 

Moreover, we have the estimates 
2N—1 n 2N 


_ 2 poo 

sup d{E(t) + y d{z(t) 

P'0 t>0 4N—2j—l f^J 0 


4N-2j 


dt 


2 N 


< 

r^j 


4jV— 2j—l 


(F.16) 


|S(0)||L-1+£K S (°) 

3 = 1 

Finally, we record a similar extension for use on q. The proof follows from an argument like 
that used in Proposition IF.21 and is thus omitted. 

Proposition F.4. Suppose that p(0) € H iN (Q) and that for j = 1,...,2 N we have that 
&lp{ 0) € H 4N ~ 2 i +1 (n). Then there exists an extension p, achieving the initial data, such that 

p € L°°([0, oo); H 4N (Q)) D L 2 ([0, oo); H 4N+1 / 2 (Q)), 
d tP € T°°([0, oo); h an ~ 1 (fi)) fl T 2 ([0, oo); H 41 ^ -1 ^ 2 (Q)), 
dipeL^([0,ocfH AN - 2 ^\n))nL 2 ([0,ooy,H 4N ~ 2 ^ +2 (n)) for j = 2,..., 2 N, 
d 2N+1 peL 2 ([0,ooy,H°m 
Moreover, we have the estimate 


sup 
t> o 


2 N 

In+ y\\ d tp( t ) 

3 = 1 


4Y-2j+l 


4AT-2j+2 


2AT+1 

+ / I lb(t)|lLv+l /2 + II^P( i )ll4iV-l/2 + y | dtP{t) 

3 =2 

2N 

SlbMIlL + Elk'rto) 

j=l 

Acknowledgements 


dt 


4N-2j+l 


(F.18) 
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